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PROPERTIES  OF  UREA-DOPED  ICE 
IN  THE  CRREL  TEST  BASIN 


Ken-ichi  Hirayama 


INTRODUCTION 

Design  of  hydraulic  structures  to  be  erected  in  re¬ 
gions  with  severe  winters  must  take  into  considera¬ 
tion  potentially  large  ice  forces.  Examination  of  past 
records  of  measured  ice  thrusts  on  similar  structures, 
and  structural  failure  due  to  ice,  is  one  way  to  esti¬ 
mate  the  ice  forces  to  be  expected.  However,  model 
studies  are  often  required  when  new  designs  are  being 
considered  or  when  the  influence  of  various  param¬ 
eters  (structure  geometry,  ice  characteristics,  etc.) 
on  ice-structure  interaction  is  being  investigated 
systematically.  Such  model  studies  and  their  quali¬ 
tative  as  well  as  quantitative  results  can  lead  to  better 
understanding  of  ice  effects  on  structures.  It  is  im¬ 
portant  that  the  model  studies  be  complemented  by 
field  measurements  and  observation  once  a  particular 
structure  has  been  built,  to  verify  model  predictions 
and  improve  modeling  techniques. 

Simple  dimensional  analysis  shows  that  physical 
modelingof  ice-structure  interaction  where  ice  fails 
in  bending  requires  that  both  the  flexural  strength 
o  and  the  elastic  modulus  E  of  ice  be  scaled  in  the 
same  ratio  as  the  geometric  characteristic  length. 
Therefore  the  ratio  El  a  must  be  the  same  for  the 
model  ice  and  the  prototype  ice.  However,  mathe¬ 
matical  analysis  of  ice  forces  on  two-dimensional 
structures  has  shown  that,  under  certain  circumstances 
which  often  prevail  in  nature,  the  peak  force  exerted 
by  the  ice  is  independent  of  E  (see,  e.g..  Yean  et  al. 
1981),  and  formulas  derived  by  Ralston  (1977), 
Danysand  Bercha  (1976)  and  Nevel  (1977)  do  not 
include  E  directly.  On  the  other  hand,  that  part  of 
the  total  resistance  of  an  icebreaker  which  is  due  to 
actually  breaking  the  ice  has  been  shown  to  be,  in 
its  dimensionless  form,  proportional  to  £/a  (Enkvist 


1972).  Some  moc  '  ais  attempt  to  correct  model 
results  for  differences  in  El  a  between  model  ice  and 
real  ice  when  the  former  falls  below  about  1500. 
Values  of  model  ice  E/o  larger  than  1500  ensure  that 
the  size  of  the  broken  floes  is  adequately  modeled. 
The  size  of  broken  floes  is  directly  related  to  8,  the 
ice  characteristic  length,  which  is  proportional  to 
the  fourth  root  of  E,  so  that  a  100%  error  in  E  leads 
to  only  a  20%  error  in  8.  Another  argument  which 
may  be  advanced  for  relaxing  the  theoretical  require¬ 
ment  of  keeping  the  ratio  E/a  the  same  for  the  mod¬ 
el  and  the  prototype  is  that  the  icebreaking  compon¬ 
ent  of  the  total  resistance  is  usually  relatively  small. 
Nevertheless,  values  of  E  for  the  model  ice  should 
remain  sufficiently  high  to  avoid  excessive  plastic 
deformation  of  the  ice  and  consequently  unrealis¬ 
tically  high  ice  forces  on  structures. 

Test  basins  without  refrigeration  capability  use 
artificial  model  ice  manufactured  from  various  mixes 
of  wax,  plastic  beads  and  other  materials.  Except 
for  that  proposed  by  Tryde  (1977),  the  exact  com¬ 
position  of  most  artificial  ice  is  proprietary  and  not 
readily  available.  Until  recently,  refrigerated  test 
basins  grew  their  model  ice  from  a  saline  solution. 
However,  this  has  restricted  model  testing  to  scales 
no  smaller  than  1  /20  because  low-strength  saline  ice 
becomes  plastic.  Furthermore,  saline  ice  leads  to 
serious  corrosion  problems  in  equipment  and  instru¬ 
mentation,  and  correspondingly  expensive  mainten¬ 
ance  and/or  protective  measures.  Timco(1979,  1980, 
1981)  proposed  that  model  ice  be  grown  from  an 
aqueous  solution  of  urea.  His  initial  studies  showed 
that  such  model  ice  could  be  used  in  studies  with  a 
geometric  scale  ratio  down  to  1/40.  One  advantage 
of  urea-doped  ice  is  the  near  elimination  of  corro¬ 
sion  problems  in  the  test  facility. 


Urea-doped  ice  has  been  used  as  model  ice  in  the 
test  basin  at  die  Ice  Engineering  Facility  of  the  U.S. 
Aimy  Cold  Regions  Research  and  Engineering  Lab¬ 
oratory  (C'RREL)  since  1980.  This  report  presents 
the  results  of  a  systematic  study  of  the  properties 
of  urea-doped  ice  based  on  measurements  on  about 
50  ice  sheets.  The  influence  of  urea  concentration, 
ice  thickness,  and  warm-up  time  (as  explained  in  the 
following  section)  on  the  mechanical  properties  of 
the  model  ice  was  investigated.  In  particular,  four 
urea  concentrations  in  the  test  basin  were  studied. 

Prior  to  August  1981  the  urea  concentration  by 
weight  in  the  test  basin  water  was  1 .05%.  In  Aug¬ 
ust  1981  the  basin  was  emptied  for  clean-up  and 
maintenance,  after  which  three  urea  concentrations 
were  successively  used,  namely  0.45%,  0.70%  and 
0.95%.  The  0.95%  concentration  was  selected  for 
permanent  use  in  the  basin.  The  study  resulted  in 
a  consistent  method  of  growth  of  an  ice  sheet  to 
achieve  any  required  ice  thickness  from  about  2  to 
8  cm,  with  bending  strength  as  low  as  about  1 5  kPa. 
Besides  ensuring  reliable  and  reproducible  ice  char¬ 
acteristics,  this  method  allows  better  scheduling  of 
test  programs.  The  particulars  of  the  method  de¬ 
scribed  in  this  report  are  valid  only  for  the  CRREL 
test  basin,  since  they  depend  upon  the  specific  re¬ 
frigeration  characteristics  of  the  facility  and  the 
urea  concentration  used.  However,  similar  methods 
can  be  devised  for  any  ice  testing  facility. 

EXPERIMENTAL  FACILITY  AND  PROCEDURES 
Ice  test  basin 

The  CRREL  test  basin  has  a  useful  length  of  34.4 
m  (1 13  ft),  a  width  of  9  m  (30  ft),  and  a  water  depth 
of  2.4  m  (8  ft).  A  5.4  m  (18  ft)  long,  3  m  (10  ft)  wide, 
2.4  m  (8  ft)  deep  preparation  tank  is  located  at  the 
east  end  of  the  main  basin.  Melting  and  recircula  ■ 
tion  tanks  are  located  at  the  west  end.  The  whole 
test  basin  is  housed  in  an  insulated,  refrigerated 
room.  The  prep-tank  area  can  be  isolated  from  the 
main  basin  by  a  vertically  sliding  insulated  door. 

The  ceiling  of  the  refrigerated  room  is  7.2  m  (24  ft) 
above  the  maximum  water  level  in  the  tank. 

Room  refrigeration  is  provided  by  seven  forced 
air  heat  exchangers  suspended  from  the  ceiling  at  a 
height  of  5.2  m  (1 7  ft)  above  maximum  water  level; 
the  total  refrigeration  capacity  is  80  kW  (2.8x  10s 
Btu/hr),  and  the  minimum  air  temperature  which 
can  be  achieved  is  approximately  -23°C  (-10°F). 

The  refrigeration  fluid  is  ammonia. 

Ice  growth  procedure 

Ice  sheets  are  grown  on  the  tank  surface  from  an 


aqueous  solution  of  urea.  The  water  is  cooled  pri¬ 
marily  by  heat  exchange  at  the  air/water  interface. 
However,  for  water  temperature  above  approximately 
2°C  (35° F),  additional  cooling  is  provided  by  pipes 
carrying  ammonia  located  in  a  trench  along  the  east 
end  of  the  tank.  At  water  temperature  below  2°C. 
these  cooling  pipes  become  fully  encased  in  ice  and 
are  ineffective.  Uniformity  of  water  temperature  in 
the  tank,  essential  for  achieving  an  ice  sheet  of  uni¬ 
form  thickness  and  mechanical  properties,  is  ensured 
by  continuously  circulating  the  water  longitudinally 
with  pumps  and  by  vertically  mixing  it  with  air  bub¬ 
bler  lines  laid  on  the  bottom  of  the  tank.  During 
the  study,  the  water  temperature  was  monitored 
regularly  at  six  locations  along  the  tank  perimeter 
by  a  thermistor  with  a  1/50°C  resolution. 

When  the  water  temperature  is  within  0.2°C  of 
the  freezing  point  of  the  urea  solution,  the  circula¬ 
tion  pumps  and  air  bubblers  are  deactivated.  Once 
the  water  surface  is  sufficiently  calm,  the  ice  sheet 
is  initiated  through  wet  seeding  by  spraying  a  fine 
mist  of  water  in  the  air,  which  is  kept  at  a  tempera¬ 
ture  of  about  -10°C  (15°F).  The  fine  mist  turns 
into  ice  crystals  which  settle  on  the  water  surface 
where  they  initiate  the  ice  sheet.  Wet  seeding  is 
necessary  to  ensure  uniform  crystal  size  in  the  ice 
sheet  (and  therefore  uniformity  of  mechanical  prop¬ 
erties)  and  a  sufficiently  high  ratio  of  ice  thickness 
to  crystal  size.  The  crystal  size  of  an  ice  sheet  grown 
in  this  manner  is  typically  of  the  order  of  2  mm  at 
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Figure  1.  Fxamples  of  measurements  of  ice  charac¬ 
teristics  ( 25-28  A  ugust  1 981  j. 


the  top  of  the  sheet  and  increases  in  the  direction 
of  ice  growth  to  about  5  mm  at  the  bottom  of  a  7- 
cm-thick  ice  sheet.  During  ice  growth,  the  air  tem¬ 
perature  in  the  room  is  kept  at  -18°  to  -10°C  (0° 
to  15°F). 

Ice  sheets  grown  at  the  above  air  temperatures 
usually  have  too  high  a  strength  for  tests  at  a  model 
scale  of  1/10  and  lower.  When  the  ice  thickness  is 
near  the  required  value  for  a  particular  test,  the  ice 
is  tempered  by  raising  the  air  temperature  in  the 
room  to  approximately  2°C  (35°F)  for  long  enough 
to  reduce  the  ice  strength  to  the  required  value.  This 
tempering  or  warm-up  technique  was  first  introduced 
by  Schwarz  (1975)  for  saline  ice.  During  the  study 
the  ice  thickness,  ice  strength,  elastic  modulus,  sur¬ 
face  temperature,  and  water  temperature  were  mon¬ 
itored  at  regular  intervals  throughout  the  warm-up 
period.  Examples  of  such  measurements  are  shown 
in  Figure  1. 


1-kg-capacity  Chatillon  push-pull  gauge  and  occasion¬ 
ally  by  a  motorized  50-lb-capacity  load  cell.  For 
about  ten  ice  sheets,  beam  tests  were  performed 
with  the  load  applied  both  upward  and  downward 
for  comparison  purposes.  For  the  other  ice  sheets, 
the  load  was  applied  only  downward. 

Elastic  or  strain  modulus 

In  the  present  study,  the  elastic  or  strain  modulus 
of  ice  E  was  measured  by  the  plate  deflection  method 
recommended  by  Sodhi  et  al.  (in  press).  Uniform 
loads  were  applied  in  incremental  steps  A P  over  a 
circular  area  of  radius  r  near  the  center  of  the  ice 
sheet,  and  the  resulting  incremental  deflection  A6 
of  the  ice  sheet  at  the  center  of  the  load  zone  was 
measured  by  an  LVDT.  The  characteristic  length  6 
of  the  ice  and  its  elastic  modulus  F.  were  calculated 
according  to  the  theory  of  an  infinite  plate  on  an 
elastic  foundation: 


Measurements 


Temperature 

As  mentioned  in  the  preceding  section,  water 
temperature  was  measured  with  a  1/50°C  resolution 
thermistor.  The  air  temperature  in  the  room  was 
continuously  monitored  and  recorded  by  two  ther¬ 
mocouples  located  immediately  below  the  heat  ex¬ 
changers.  During  ice  growth  and  tempering,  the  ice 
surface  temperature  was  measured  by  thermistors 
mounted  on  a  piece  of  Styrofoam. 

Ice  thickness 

During  ice  growth  the  ice  thickness  was  regularly 
measured  with  a  caliper  with  0.1  mm  accuracy  at 
six  points  along  the  perimeter  of  the  test  basin. 
After  a  test,  the  ice  thickness  near  the  track  left  by 
the  tested  structure  was  also  measured  to  determine 
the  distribution  of  ice  thickness  and  to  evaluate  its 
degree  of  uniformity. 

Ice  flexural  strength 

The  ice  flexural  strength  a  was  determined  from 
in-situ  cantilever  beam  tests.  The  beam  width  B 
was  taken  as  one  to  two  times  the  beam  thickness 
H,  and  its  length  L  was  6  to  7  times  the  thickness 
to  minimize  the  buoyancy  effects  (Tatinclaux  and 
Hirayama  1982)  so  that  the  formula  for  simple  can¬ 
tilever  beams  could  be  applied : 


a  =  — 
BH2 


(1) 


where  yw  is  the  specific  weight  of  water,  w  is  the 
Poisson’s  ratio  of  ice,  taken  equal  to  1  /3,  a  =  r/K 
and  lny  is  the  Euler’s  constant. 

This  method,  which  is  nondestructive  and  simple 
to  apply,  proved  to  give  consistent  results.  The  set¬ 
up  for  load  application  and  deformation  measure¬ 
ment  either  in  the  upward  or  downward  mode  is 
shown  schematically  in  Figure  2. 

Urea  concentration 

Urea  concentration  in  the  tank  water  and  in  the 
ice  was  measured  by  a  Beckman  Blood-Urea-Nitrogen 


where  P  is  the  failure  load  applied  at  the  free  end 
of  the  beam.  The  load  was  usually  applied  by  a 
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Figure  2.  Experimental  set-up  for  measurements  of 
elastic  modulus. 


(BUN)  Analyzer  2,  which  had  been  initially  cali¬ 
brated  by  urea  solutions  in  distilud  water  of  known 
concentration. 


ICE  GROWTH  AND  STRUCTURE 

Ice  thickness  distribution 

Two  examples  of  ice  thickness  distribution  in  the 
test  basin  are  shown  in  Figure  3.  In  these  two  ex¬ 
amples,  taken  during  the  summer  months,  besides 
local  effects  from  the  melting-recirculation  tanks 
at  the  north  end  of  the  tank  and  from  the  prep  room 
at  the  south  end,  the  ice  thickness  can  be  observed 
to  increase  slightly  from  the  northwest  corner  to 
the  southeast  corner  of  the  basin.  This  trend  was 
reversed  during  the  winter  season;  that  is,  the  ice 
increased  in  thickness  from  the  southeast  corner  to 
the  northwest  comer  of  the  basin.  The  outsides  of 
the  north  and  west  walls  are  exposed  to  the  atmos¬ 
phere,  while  the  south  wall  separates  the  basin  from 
the  instrumentation  corridor,  where  air  conditioning 
keeps  the  air  temperature  at  20°-25°C  year-round. 


Reasons  for  the  slight  nonuniformity  in  the  ice  thick¬ 
ness  which  have  been  identified  or  are  suspected  are: 

1 .  Uneven  air  temperatures  outside  the  basin 
room  (atmosphere,  instrumentation  cor¬ 
ridor,  prep  room). 

2.  Uneven  temperature  distribution  and  nat¬ 
ural  convection  in  the  basin  room  due  to 
inadequate  or  uneven  defrosting  of  the 
heat  exchanger  coils. 

3.  Excessive  traffic  in  the  room,  resulting  in 
excessive  opening  and  closing  of  the  doors. 

4.  Excessive  heating  of  meltwater  in  the 
melting  tank. 

5.  Inadequate  operation  of  the  pump  recir¬ 
culation  system. 

6.  Nonuniform  air  bubbling  prior  to  seeding. 
Reasons  2-6  can  be  greatly  reduced  if  not  totally 
eliminated  by  careful  operation  of  the  ice  test  basin. 
In  addition,  it  was  observed  that  the  longitudinal 
center  portion  (4  to  32  m)  of  the  basin,  where  the 
model  tests  are  usually  carried  out,  has  a  satisfactory 
uniform  ice  thickness  distribution.  Table  1  presents 
statistics  of  ice  thickness  distribution  in  the  test 


Figure  3.  Examples  of  ice  thickness  distribution  in  tank. 


Table  1 .  Statistical  characteristics  of  ice  thickness  distribution. 


Ice  thickness  Coef.  of 


Date 

Measurement 

zone 

No.  of 
points 

Mean 

(cm) 

Std.  dev. 
(cm) 

max. 

(cm) 

min. 

(cm) 

variation 

(%) 

28  Aug 

6—30  m 

27 

4.14 

0.176 

4.40 

3.77 

4.3 

6  July 

6-30  m 

27 

7.29 

0.172 

7.58 

6.95 

2.4 

9  Sept 

6-27  m 

24 

6.4  S 

0.160 

6.65 

6.04 

2.5 

14  Sept 

6-27  m 

24 

3.42 

0.047 

3.55 

3.37 

1.4 

IS  Sept 

6-24  m 

21 

3.37 

0.083 

3.57 

3.20 

2.5 

24  Sept 

8-28  m 

20 

4.67 

0.070 

4.80 

4.55 

1.5 

16  Nov 

8-30  m 

12 

$.61 

0.090 

5.74 

5.42 

1.6 

t9  Nov 

8-30  m 

11 

3.29 

0.053 

3.39 

3.22 

1.6 

23  Nov 

8-30  m 

12 

5.61 

0.047 

5.67 

5.53 

0.8 

27  Nov 

8-30  m 

12 

3.41 

0.078 

3.55 

3.26 

2.3 

1  Dec 

8-30  m 

11 

5.55 

0.033 

5.58 

5.42 

0.6 

4  Dec 

1 0-30  m 

11 

3.35 

0.050 

3.47 

3.27 

1.5 

1.9% 


h,  Average  of  Ice  Thickness  at  6  Points  (cm) 

Figure  4.  Average  ice  thickness  in  test  area  of  tank 
versus  average  thickness  at  six  points  along  the 
tank  boundary. 

portion  of  the  basin.  For  most  ice  sheets  the  stan¬ 
dard  deviation  in  ice  thickness  is  less  than  1  mm, 
and  of  the  order  of  2%  or  less  of  the  mean  thickness. 
Finally,  as  shown  in  Figure  4,  the  mean  ice  thick¬ 
ness  hQ  in  the  test  portion  of  the  basin  was  practi¬ 
cally  equal  to  the  average,  h,  of  the  thicknesses 
measured  at  the  six  points  at  the  periphery  of  the 
basin  shown  in  the  insert  of  Figure  1 .  In  the  re¬ 
mainder  of  the  report,  this  average  thickness  has 
been  used  as  the  nominal  thickness  of  the  ice  sheet. 

Ice  growth  during  freeze-up 

During  ice  growth  the  water  at  the  tank  surface 
and  the  air  in  the  room  are  practically  quiescent. 


Therefore  the  classical  relationship  between  ice 
thickness  and  negative  degree-hours,  found,  for  ex¬ 
ample,  in  Ashton’s  (1978)  or  Calkins’s  (1979)  pa¬ 
pers,  for  lake  ice  should  be  applicable,  namely 


(4) 


where 

kt  -  heat  conductivity  of  ice 
hia  =  heat  transfer  coefficient  at  ice/air  inter¬ 
face 

P;  =  ice  density 
X  =  latent  heat  of  fusion  of  ice 

E=P  (Ta-Tiw)dt 

o 


(Ta  =  air  temperature,  Tiw  =  freezing  temperature 
of  urea  solution,  t  =  time).  Equation  4  assumes  that 
there  is  no  temperature  gradient  in  the  meltwater 
and  that  the  temperature  distribution  in  the  ice  sheet 
is  linear,  as  shown  in  Figure  5.  In  Figure  6,  the 
thicknesses  of  the  ice  sheets  studied  are  plotted 
versus  (-2)  in  °C-hr,  where 

E=f  Tadt. 

o 

That  is,  Tiw,  the  ice/water  interface  temperature  or 
freezing  temperature  of  the  urea  solution,  was  taken 
equal  to  0°C,  neglecting  depression  of  the  freezing 
point  of  the  urea  solution  (0.1 55°C  for  a  0.5%  solu¬ 
tion  and  0.3 1°C  for  a  1%  solution)  with  respect  to 
the  freezing  point  of  water.  The  data  were  separated 
into  two  sets.  The  first  set  (Fig.  6a)  was  obtained 
prior  to  August  1981  and  the  second  (Fig.  6b)  after 
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Figure  5.  Definition  sketch  for  heat  transfer  through  a 
floating  ice  sheet. 
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Figure  6.  Relationship  between  ice  thickness  and 
degree-hours. 


August  1981 ,  when  control  of  water  temperature  at 
seeding  and  frequency  of  heat  exchanger  defrosting 
were  improved.  As  can  be  seen  from  the  two  fig¬ 
ures,  these  improvements  resulted  in  a  slightly  larger 
growth  rate  of  the  ice  and  in  a  significant  reduction 
in  the  scatter  of  the  data  points.  Lines  correspond¬ 
ing  to  equation  4  were  fitted  through  the  data  points 
in  Figure  6.  In  Figure  6a,  the  heat  conductivity  of 
ice  was  assumed  to  be  that  of  freshwater  ice,  i.e. 
k{  =  19.26  cal/°C-cm-hr,  and  hia  was  calculated  by 
nonlinear  regression  analysis  of  the  experimental 
data  and  found  equal  to 

hia  =  0.73  cal/hr-cm2-°C  . 

In  Figure  6b,  both  ki  and  hia  were  calculated  by 
nonlinear  regression  analysis  to  be 

fcj  =  14.85  cal/°C-cm-hr 

hia  =  0.81  cal/°C-cm2-hr . 

In  these  calculations  it  was  assumed  that  X  =  79.8 
cal/g  and  p{  =  0.92  g/cm3. 

As  also  shown  in  Figure  6,  for  ice  thickness  less 
than  2  to  3  cm,  eq  4  can  be  approximated  by  a  lin¬ 
ear  relationship 

*"j&E  ® 

as  proposed  by  Assur  and  Weeks  (1964). 

Ice  growth  during  warm-up 

It  was  observed  that  during  the  warm-up  or  tem¬ 
pering  period  necessary  to  achieve  the  ice  strength 
required  for  a  particular  model  scale,  the  ice  thickness 


Figure  7.  Relative  increase  of  ice 
thickness  during  warm-up. 


continued  to  increase  slightly.  Let  hx  be  the  ice 
thickness  at  the  start  of  the  warm-up  period  when 
room  refrigeration  is  discontinued  and  the  air  tem¬ 
perature  is  allowed  to  rise  above  freezing,  and  h(t) 
the  ice  thickness  at  any  subsequent  time.  The  rela¬ 
tive  increase  in  thickness 

M  _  h{t)-h  j 
h;  =  h; 

has  been  plotted  versus  t/hf  in  Figure  7.  The  param¬ 
eter  r//ij2  was  selected  on  the  basis  of  the  study  by 
Schwarz  and  Miloh  (1972),  who  showed  that  the 
time  necessary  for  a  saline  ice  sheet  initially  at  uni¬ 
form  temperature  to  reach  equilibrium  temperature 
when  one  face  of  the  ice  is  suddenly  brought  to  a 
different  temperature  TA  is  independent  of  TA  and 
of  the  ice  salinity  but  inversely  proportional  to  the 
square  of  the  ice  thickness.  In  spite  of  the  scatter 
in  the  data  points.  Figure  7  shows  that  ^h/hi  in¬ 
creases  asymptotically  toward  a  value  of  approxi¬ 
mately  0.20.  The  equation  of  the  curve  in  Figure  7 
was  obtained  by  nonlinear  regression  analysis  as 

M  0.22  tfh] 

hx  t/hf  +  0.57  ' 

The  final  ice  thickness  hf  at  which  the  model  tests 
were  conducted  is  plotted  versus  the  initial  ice  thick¬ 
ness  /ij  in  Figure  8.  It  can  be  seen  that  the  data  points 
lie  on  a  straight  line: 

hf  =  1.17  h{ .  (7) 


Figure  8.  Final  ice  thickness  versus  initial  ice  thick¬ 
ness  at  start  of  warm-up. 
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Because  there  is  less  scatter  in  Figure  8  than  in  Fig¬ 
ure  7,  eq  7  may  be  more  reliable  than  eq  6  for  pre¬ 
dicting  the  ice  thickness  at  which  warm-up  should 
be  started  to  achieve  a  desired  final  ice  thickness  as 
long  as  the  test  conditions  are  within  the  limits  of 
those  investigated  in  the  present  study. 

Structure  of  urea-doped  ice 

Photographs  of  thin  sections  of  urea-doped  ice 
are  shown  in  Figure  9.  It  can  be  seen  that  urea- 
doped  ice  is  similar  to  natural  sea  ice  insofar  as  it 
contains  distributed  brine-like  cells  (see  Fig.  9b). 
However,  contrary  to  sea  ice,  the  urea-doped  ice 
grown  in  the  CRREL  basin  is  composed  of  two  dis¬ 


tinct  layers,  as  evidenced  in  Figure  9a.  The  thicker 
bottom  layer  has  a  c-axis-horizontal,  columnar  struc¬ 
ture  characteristic  of  dendritic  growth,  while  the 
thinner  top  layer  is  made  of  much  shorter  columns 
with  mixed  vertical  and  horizontal  c-axis  orientation. 
As  is  obvious  in  Figure  9a,  there  is  a  distinct  discon¬ 
tinuity  in  crystal  structure  at  the  interface  between 
the  top  and  bottom  layers.  The  thickness  of  the 
top  layer  of  urea-doped  ice  is  relatively  much  larger 
than  that  of  the  transition  incubation  layer  in  normal 
sea  ice. 

The  growth  of  urea-doped  ice  may  be  compared 
to  the  congelation  process  of  metal  alloy  melts  as 
studied  by  Tiller  (1958).  If  the  analogy  is  valid. 


according  to  Tiller  columnar  growth  of  urea  ice 
takes  place  when  the  following  instability  condition 
is  satisfied: 


D  .  T  GDk  1 
kR  "  ‘  c0Rm(\  -*)J  '  (  ’ 


mc0  (1  ~k ) 
Dk 


G  =  temperature  gradient  in  the  liquid  at  the 
freezing  interface 
R  =  growth  rate 

k  =  equilibrium  partition  coefficient  of  the 
solute 

m  =  liquidus  slope  (a  negative  quantity) 

D  -  diffusion  coefficient 
c0  =  initial  concentration  of  solution. 

The  thickness  hx  of  the  top  layer  with  random  c-axis 
orientation  of  the  crystals  is  given  by  Tiller  (1 958) 
as 


It  is  uncertain  whether  this  equation  is  applicable 
to  the  case  of  very  thin  urea-doped  ice  where  con¬ 
nections  between  brine  cells  are  possible  and  brine 
migration  is  to  be  expected.  Moreover,  it  has  been 
observed,  as  discussed  in  the  following  section,  that 
the  thickness  of  the  top  layer  of  model  ice  continues 
to  increase  during  the  growth  of  the  ice  sheet.  This 
puzzling  phenomenon  remains  unexplained,  but 
would  imply  a  recrystallization  process  at  the  inter¬ 
face  between  the  top  and  bottom  layers  of  the  model 
ice.  The  increase  in  top  layer  thickness  with  increas¬ 
ing  total  ice  thickness  is  illustrated  in  Figure  10a  for 
the  1.05'/f  urea  concentration  and  in  Figure  I  Ob  for 
the  other  three  concentrations  investigated.  In  these 
figures  it  is  evident  also  that  the  ratio  hJH  decreases 
with  increasing  urea  concentration,  and  decrc?ses 


Ttmporotur*  of  Wotfr  of  S4M0i*9  (^C) 

Figure  II.  Effect  of  water  temperature  at  time  of  seeding  on 
top  layer  thickness  (c  =  !.057o). 

Table  2.  Urea  concentration  in  ice  and  water  solution. 


Date  Time 


tee  thickness  Concentration  (%) 

h  (cm)  Top  Bottom  Solution  Remarks 


20  Nov 

21  Nov 


27  Nov 

28  Nov 

29  Nov 

30  Nov 


Warm-up 


Warm-up 


with  increasing  thickness  h  toward  an  asymptotic 
value  varying  between  0.2  and  0.3.  One  other  param¬ 
eter  which  was  found  to  affect  the  top  layer  thick¬ 
ness  is  the  water  temperature  at  seeding.  Figure  1 1 , 
for  the  particular  urea  concentration  of  1 .05%,  shows 
an  increase  in  the  ratio  hjh  with  increasing  water 
temperature  at  seeding.  The  observed  variations  of 
h(  with  urea  concentration  and  water  temperature 
at  seeding  are  in  qualitative  agreement  with  eq  9. 

Finally,  because  of  the  difference  in  their  crystal 
structure,  the  urea  concentrations  in  the  top  and 
bottom  layers  are  expected  to  be  different.  This  was 
confirmed  by  measurements  listed  in  Table  2,  where 
it  can  be  seen  that  during  the  freeze-up  period  the 
urea  concentration  in  the  top  layer  is  of  the  order  of 
half  that  in  the  bottom  layer,  itself  about  half  the 
urea  concentration  in  the  basin  water.  Once  warm¬ 
up  of  the  ice  sheet  was  started,  the  urea  concentration 
in  both  layers  was  observed  to  decrease,  which  indi¬ 
cates  an  increase  in  cavities  within  the  ice,  leading  to 
drainage  and  loss  of  urea. 


MECHANICAL  PROPERTIES  OF 
UREA-DOPED  ICE 

Introductory  remarks 

Studies  on  the  mechanical  properties  of  saline 
ice,  namely  bending  strength  and  strain  modulus, 
have  shown  that  these  properties  are  functions  of 
many  parameters,  such  as  rate  of  loading,  solute  con¬ 
centration  in  the  meltwater,  ice  temperature,  direc¬ 
tion  of  load  application,  and  even  method  of  meas¬ 
urement.  In  the  present  study,  all  the  measurements 
were  made  in  situ.  As  mentioned  earlier,  the  bending 
strength  was  measured  in  small  cantilever  beam  tests, 
and  the  strain  modulus  by  the  plate  deflection  meth¬ 
od.  In  the  majority  of  the  cantilever  beam  tests,  the 
load  was  applied  manually  and  at  such  a  rate  that 
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Figure  12.  Relationship  between  ice  sur¬ 
face  temperature  and  air  temperature. 
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beam  failure  occurred  within  1  to  2  seconds,  follow¬ 
ing  the  recommendations  of  the  1AHR  Committee 
on  Ice  Problems.  Since  the  tests  were  performed  in 
situ,  the  bottom  of  the  ice  was  always  at  melting 
temperature,  and  only  the  ice  surface  temperature 
could  be  varied,  by  varying  the  ambient  air  tempera¬ 
ture.  Figure  12  shows  that  a  variation  of  some  I2°C 
(22°F)  in  air  temperature  results  in  a  variation  of 
only  about  2°C  (3.6°F)  in  the  ice  surface  temperature. 

In  addition  to  the  effects  of  air  temperature,  urea 
concentration  and  direction  of  loading  on  ice  proper¬ 
ties,  the  influence  of  the  two-layer  structure  of  the 
ice  on  the  model  ice  properties  can  be,  at  least  qual¬ 
itatively,  predicted  by  a  mathematical  model,  pre¬ 
sented  below. 

Last  but  not  least  is  the  effect  of  tempering  or 
warm-up  on  the  mechanical  properties  of  urea-doped 
ice.  The  effect  of  warm-up  was  investigated  at  great 
length  in  order  to  be  able  to  predict  the  warm-up 
duration  necessary  to  achieve  the  required  mechan¬ 
ical  properties,  primarily  bending  strength,  and  con¬ 
sequently  to  establish  in  advance  a  reliable  schedule 
for  ice  freeze-up,  ice  tempering,  and  ultimately  the 
test  program. 

The  experimental  data  on  the  mechanical  proper¬ 
ties  of  untempered  urea-doped  ice  have  been  tabu¬ 
lated  in  Appendix  B  and  those  for  tempered  ice  in 
Appendix  C. 

Model  of  a  two-layer  elastic  material 

Consider  a  sheet  of  thickness  H  of  a  nonhomo- 
geneous  material  composed  of  two  layers.  The  top 
layer,  of  thickness  aH ,  is  a  homogeneous,  elastic 
material  characterized  by  its  bending  strength  a, 
and  elastic  modulus  E , .  The  bottom  layer,  of  thick¬ 
ness  ( !  - a)H ,  is  characterized  by  o2  and  E2 ,  with 
o,  >  a2  and  E ,  >  E2  (see  Fig.  13). 

A  cantilever  beam  of  length  A  and  width  B  of 
this  composite  material  is  subjected  to  a  load  P  at 
the  free  end.  It  is  assumed  that  the  strain  distribu¬ 
tion  e(z)  across  the  beam  is  linear,  as  shown  in  Figure 
13,  i.e. 

<*z)=e,z/z0  (10) 

where  e,  is  the  strain  at  the  top  surface  of  the  beam, 
z  is  measured  positive  upward  from  the  yet-unde- 
termined  neutral  axis,  and  z0  is  the  distance  from 
the  neutral  axis  to  the  beam  top  surface. 

The  stress  distribution  in  the  top  layer  is  given  by 

0*1)  =  ^6  =  £',€,  7- 
"0 
or 

o(l)  =  o,  7  (r0  ~ aff  <  z  <  z0)  (II) 

zo 


oh 

(l-a)H  Eg ,  0% 

p 


Width  B 
Thickness  H 
Length  L 


Figure  13.  Definition  sketches  for  two-layer  model. 


where  a,  is  the  stress  at  the  top  surface.  In  the  bot¬ 
tom  layer  the  stress  distribution  is  given  by 

E.  z 

o<2)  =  E2e  =  —  a,  —  (z0-aH  <  z  <  H). 
bi  zo 

(12) 


From  the  condition  /  adz  =  0,  the  position  of  the 
neutral  axis  is  found  to  be 


H  1  - a 2  +a2b 
Z°  2  \  -a  +ab 


(13) 


where  b  =  jE2 . 

If  E0  is  the  elastic  modulus  of  the  equivalent  homo¬ 
geneous  beam,  i.e.  which  experiences  the  same  tip 
displacement  6  under  the  same  load  P ,  then 


EqI  / j  +  E j  12 


(14) 


where  /  =  BH31 12  is  the  beam  moment  of  inertia 
about  its  mid-thickness  and  /,  and  /2  are  the  mo¬ 
ments  of  inertia  of  the  top  layer  and  bottom  layer, 
respectively,  about  the  actual  neutral  axis.  Equation 
14  yields 


f®  =aj+  (L-fl)3  +  jgliigi 

£,  b  1  -a  +  ob  ' 


(15) 


Application  of  eq  IS  to  measurements  carried  out  on 
urea-doped  ice  beams  and  their  top  layers  alone  (i.e. 
after  scraping  off  the  bottom  layer)  yielded  values  of 
£■,  /E2  of  20  or  more. 

When  the  load  P  is  applied  downward,  failure 
occurs  when  the  stress  at  the  top  of  the  beam  o,  be¬ 
comes  equal  to  the  tensile  strength  o,  of  the  top 
layer  material.  Let  o0  be  the  nominal  strength  of 
the  equivalent  homogeneous  beam  with  neutral  axis 
at  mid-thickness,  i.e.  calculated  from  eq  I  with  h=  H, 
it  can  be  shown  that 


(16) 


When  the  load  P  is  applied  upward,  it  can  be  ex¬ 
pected  that  a  crack  forming  at  the  bottom  of  the 
beam  when  the  stress  there  becomes  equal  to  a2  will 
propagate  up  to  the  interface  between  the  top  and 
bottom  layers.  Total  failure  of  the  beam  will  then 
occur  when  the  stress  at  the  interface  becomes  equal 
to  o, .  The  corresponding  nominal  strength  o'0  of 
the  composite  beam  calculated  from  eq  1  with  h  =  H 
is  such  that 


Therefore,  when  such  a  two-layer  beam  is  tested  in 
both  the  downward  and  upward  modes,  the  ratio  of 
the  corresponding  strengths  given  by  application  of 
simple  beam  theory  (eq  1 )  is  given  by 


0l=2a2Z_±tjL 

On  H  En 


(18) 


The  reasoning  which  led  to  eq  17  and  18  should  be 
valid  for  urea-doped  ice,  where  the  bottom  layer,  with 
its  columnar  structure  and  the  brine  pockets  it  con¬ 
tains,  offers  little  resistance  to  tension,  especially 
during  warm-up. 

The  nominal  ratio  E0/o0  can  be  obtained  in  terms 
of  the  ratio  E ,/o,  for  the  top  layer  from  eq  13  and 
16  as 


fo/q°  =  ,  fl  =  1  -a2  +a2b 

El/ol  H  1  -a  +ab  v  ' 

Equations  13,  15,  16  and  18  have  been  plotted  in 
Figure  14  with  a  as  the  variable  and  for  the  particular 
value  Et  IE2  =  20.  It  can  be  noted  that,  according  to 
the  above  model,  over  a  fairly  wide  range  of  the  param¬ 
eters  (approximately  0.1  <  a  <  0.4)  both  the  ratios 
a0/o,  and  E0/El  remain  nearly  constant,  while  the 
ratio  o'0/o0  varies  significantly  from  0.03.  at  a  =  0.1 . 
to  0.47,  at  a-  0.4,  for  the  particular  value  of 
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Figure  14.  Variation  of  properties  of  a  two-layer  ma¬ 
terial  with  relative  thickness  of  top  layer. 


Ex  IE2  =  20.  This  range  of  a  does  include  the  values 
of  ht/h  of  0.2  to  0.3  observed  for  the  urea-doped  ice 
investigated  in  the  present  study. 

Finally,  it  can  be  inferred  from  Figure  14  that, 
according  to  the  double  layer  model,  when  a  low 
value  of  o0,  i.e.  of  o0/a, ,  is  required  together  with 
a  high  value  of  E0/o0,  i.e.  of  ( E0lo0)/(E}  /a, ),  the 
top  layer  should  be  very  thin,  of  the  order  of  5%  or 
less  of  the  total  thickness. 

Properties  of  urea-doped  ice  during  freeze-up 

The  measurements  of  ice  thickness,  and  ice  me¬ 
chanical  properties  performed  during  freeze-up,  have 
been  regrouped  in  Table  3  for  easier  reference. 

Flexural  strength 

Figure  15  presents  results  for  measurements  of 
the  flexural  strength  a  (in  the  downward  mode)  for 
ice  grown  from  a  1 .05%  urea  solution.  The  results 
pertain  to  several  ice  sheets  grown  at  various  air  tem¬ 
peratures,  and  measurements  were  taken  at  different 
ice  thicknesses  during  the  freezing  period.  In  spite 
of  the  unavoidable  scatter  in  the  data,  it  is  apparent 
that  the  flexural  strength  increases  with  decreasing 
air  temperature,  i.e.  with  decreasing  ice  surface  tem¬ 
perature.  As  mentioned  previously,  a  drop  of  about 
1 1°C  (20°F)in  air  temperature  ft  resulted  in  a  drop 


Figure  15.  Variation  of  flexural  strength  with 
ice  thickness  and  air  temperature  during  freeze- 
up  (c  =  1.05%). 


of  about  2°C  (4°F)  in  ice  surface  temperature.  This 
in  turn  translated  into  an  increase  in  ice  strength  of 
the  order  of  30  kPa,  from  50-60  kPa  at  T  >  7°C 
(20°F)  to  80-100  kPa  at  Ta  <  -12°C(10'T)foran 
ice  thickness  of  about  3  cm,  as  can  be  seen  from  Fig¬ 
ure  15. 

Figure  1 5  also  shows  an  increase  in  flexural  strength 
with  ice  thickness,  which  parallels  the  increase  in  top 
layer  thickness  with  total  thickness  (Fig.  10),  thereby 
further  indicating  the  important  effect  of  the  top  layer 
on  the  nominal  flexural  strength  of  the  ice. 

A  similar  increase  of  flexural  strength  with  ice 
thickness  was  observed  with  the  ice  sheets  grown  from 
solutions  of  varying  urea  concentration  as  evidenced 
in  Figure  16,  where  the  ice  sheets  were  grown  at  an 
air  temperature  between  -12°  and  -9°C  (10°  and 
15°F.  Average  empirical  relationships  between  ice 
strength  and  ice  thickness  were  calculated  by  a 


Figure  16.  Variation  of  flexural  strength  with  ice 
thickness  during  freeze-up  (all  urea  concentra¬ 
tions). 
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Table  3.  Summary  of  measurements  during  freeze-up 


E/o 

C%  Date  Time  o  (kPa)  a'(kPa)  E  (MPa)  E'(MPa)  h  (cm )  h  t(cm/  o'  /a  XI CT h,/h 


28  May 

0900 

44.0 

2.34 

0.72 

0.31 

1400 

34.7 

2.93 

0.88 

0.30 

28  May 

2100 

32.1 

3.16 

0.91 

0.29 

29  May 

0800 

23.0 

3.40 

1.05 

0.31 

1400 

20.3 

31.6 

1 .56 

0.31 

2  June 

0900 

67.6 

1.69 

0.61 

0.36 

1600 

83.9 

2.43 

0.84 

0.35 

3  June 

0900 

84.0 

3.72 

1.08 

0.29 

1200 

85.1 

4.01 

1.18 

0.29 

4  June 

1030 

66.6 

91.7 

1 .85 

0.68 

1.37 

0.36 

1630 

78.2 

2.38 

0.83 

0.35 

5  June 

0830 

87.1 

3.75 

1.03 

0.27 

1100 

100.6 

4.03 

1.13 

0.28 

6  June 

0800 

49.7 

1.73 

0.71 

0.41 

2030 

63.6 

2.82 

0.98 

0.35 

7  June 

0900 

66.7 

3.76 

1.10 

0.29 

1900 

67.0 

4.44 

1.25 

0.28 

8  June 

0830 

94.0 

5.34 

1.45 

0.27 

1100 

75.0 

5.61 

1.S3 

0.27 

10  June 

0900 

74.3 

2.52 

0.75 

0.30 

1630 

74.4 

62.8 

3.31 

0.87 

0.84 

0.26 

1 1  June 

0830 

90.8 

59.8 

4.43 

1.10 

0.66 

0.23 

13  June 

2300 

48.9 

43.7 

3.89 

0.88 

0.89 

0.23 

14  June 

0930 

88.1 

5.24 

1.06 

0.20 

15  June 

0930 

79.9 

7.38 

1.16 

0.16 

1200 

163.9 

16.3 

1.13 

0.16 

18  June 

0900 

53.7 

2.66 

0.72 

0.27 

1600 

60.8 

3.26 

0.91 

0.28 

19  June 

0900 

93.6 

4.82 

1.03 

0.21 

20  June 

0800 

68.2 

2.14 

0.72 

0.34 

21  June 

0800 

89.8 

4.34 

1.05 

0.24 

22  June 

0900 

90.8 

6.67 

1.31 

0.20 

25  June 

0830 

62.9 

2.21 

0.76 

0.34 

26  June 

0800 

91.9 

4.25 

1.16 

0.27 

27  June 

0800 

52.3 

1.92 

0.69 

0.36 

28  June 

0800 

81.0 

4.00 

1.00 

0.25 

29  June 

0800 

89.7 

6.01 

1.21 

0.20 

1130 

76.6 

I  July 

0800 

86.8 

3.5  7 

0.76 

0.21 

3  July 

1130 

58.0 

2.77 

0.31 

4  July 

1030 

74.8 

4.02 

1.05 

0.26 

5  July 

1330 

92.7 

5.80 

1.31 

0.23 

6  July 

0900 

123.2 

6.98 

1.45 

0.21 

12  30 

94.9 

84.6 

0.89 

0.21 

1400 

80.2 

7  4.0 

0.92 

0.21 

8  July 

1300 

90.7 

59.7 

2.91 

0.89 

0.66 

0.31 

9  July 

0800 

87.0 

93.9 

4.73 

1.10 

1.08 

0.23 

15  July 

0800 

94.0 

2.80 

1.14 

0.41 

0930 

98.6 

3.00 

1.28 

0.43 

20  July 

0900 

84.2 

3.66 

1.41 

0.39 

21  July 

0800 

88.2 

134.0 

4.71 

1.64 

1.52 

0.35 

24  July 

0830 

77.9 

5  1.0 

4.00 

1.23 

0.63 

0.31 

19  Aug 

0900 

76.2 

40.8 

2.08 

1.07 

5.35 

0.51 

25  Aug 

2030 

126.5 

85.2 

27.1 

1.60 

0.79 

0.67 

2.14 

0.49 

26  Aug 

1100 

167.1 

104.9 

33.7 

2.93 

1.22 

0.63 

2.01 

0.42 

1500 

203.6 

81.3 

34.1 

0.40 

1.67 

30  Aug 

0930 

173.4 

76.4 

220.5 

:  "6 

l.ll 

0.44 

1.27 

0.43 

1900 

193.4 

90.9 

31  1.9 

3.  / 

1.31 

0.47 

1.61 

0.37 

31  Aug 

0900 

167.1 

99.8 

306.5 

4.1 

1.60 

0.60 

1.83 

0.33 
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E  to 
i  n-3 


Date 

Time 

a  fkPa) 

o'  (k  Pa) 

E  (MPa) 

E'(MPa) 

h  (cm) 

h  x(cm) 

o'  (o 

>> 

» 

*  1 

h,/h 

6  Sept 

1000 

137.7 

74.6 

187.3 

176.5 

2.91 

1.05 

0.54 

1.36 

0.36 

1800 

133.8 

90.3 

169.3 

191.8 

3.82 

1.20 

0.68 

1.27 

0.31 

7  Sept 

0830 

171.7 

80.3 

178.8 

181.7 

5.42 

1.45 

0.47 

1.04 

0.27 

10  Sept 

2130 

102.2 

58.1 

185.8 

1.61 

0.74 

0.57 

1.82 

0.46 

1 1  Sept 

0900 

125.0 

64.5 

180.3 

189.6 

2.85 

1.04 

0.52 

1.44 

0.36 

1 5  Sept 

0930 

87.1 

50.4 

69.1 

2.02 

0.81 

0.58 

0.79 

0.40 

20  Sept 

0800 

65.1 

49.2 

68.7 

70.9 

2.76 

0.91 

0.76 

1.06 

0.33 

27  Sept 

0900 

66.5 

35.1 

79.0 

86.0 

2.38 

0.76 

0.53 

1.19 

0.32 

7  Oct 

1300 

76.4 

35.0 

78.4 

2.36 

0.76 

0.46 

1.03 

0.32 

8  Oct 

0900 

122.5 

47.6 

80.8 

4.25 

1.06 

0.39 

0.66 

0.24 

17  Oct 

2130 

81.6 

43.7 

118.4 

1.93 

0.75 

0.54 

1.45 

0.39 

18  Oct 

1600 

106.9 

53.5 

131.1 

3.21 

1.00 

0.50 

1.23 

0.31 

2000 

121.5 

56.1 

122.5 

3.75 

1.09 

0.46 

1.01 

0.29 

19  Oct 

0900 

163.5 

59.4 

158.2 

4.61 

1.32 

0.41 

1.10 

0.29 

13  Nov 

0800 

81.2 

1.99 

0.58 

0.29 

14  Nov 

1030 

124.8 

5.1  5 

1.00 

0.19 

18  Nov 

0800 

94.3 

62.2 

2.71 

0.69 

0.66 

0.2  6 

20  Nov 

1200 

87.2 

66.0 

1.99 

0.61 

0.76 

0.31 

21  Nov 

0830 

143.7 

187.9 

4.92 

1.07 

1.31 

0.22 

26  Nov 

0800 

92.1 

74.8 

2.56 

0.68 

0.81 

0.27 

28  Nov 

0900 

60.8 

45.8 

1.99 

0.73 

0.75 

0.37 

1800 

103.0 

88.7 

2.87 

0.85 

0.86 

0.30 

29  Nov 

1800 

146.2 

182.7 

4.80 

1.08 

1.25 

0.23 

2  Dec 

1200 

91.0 

64.2 

2.76 

0.73 

0.71 

0.26 

nonlinear  regression  analysis  of  the  data  presented  in 
Figure  15  (-12°C  <  Ta  <  -9°C)and  Figure  16  as 


C=  1.05% 

142  h 

0  h  +  2.84 

(20a) 

C  =  0.95% 

261  h 
°  h  +  5.06 

(20b) 

c  -  0.70% 

256  h 
°  h  +  2.88 

(20c) 

c  -  0.45% 

329/j 
°  h  +  2.61 

(20d) 

where  h  is  expressed  in  cm  and  a  in  kPa.  Curves  cor¬ 
responding  to  the  above  equations  are  drawn  in  Fig¬ 
ures  15  and  16. 

Figure  16  shows  a  significant  increase  in  flexural 
strength  with  decreasing  urea  concentration  in  the 
melt  solution.  For  example,  for  2-cm-thick  ice  the 
measured  flexural  strengths  at  urea  concentrations 
of  0.45%,  0.7%,  0.95%  and  1.05%  were  130  kPa, 

100  kPa,  70  kPa  and  50  kPa,  respectively,  and  for 
6-cm-thick  ice,  the  corresponding  flexural  strengths 
were  about  210  kPa,  150  kPa,  130  kPa  and  90  kPa, 
respectively.  To  further  illustrate  the  effect  of  urea 
concentration  c  on  bending  strength,  values  of  a 
predicted  by  eq  20  are  plotted  versus  c  for  four  ice 
thicknesses  in  Figure  1 7.  It  can  be  seen  that  a  is  a 
decreasing  linear  function  of  c.  A  nonlinear  regres¬ 
sion  analysis  of  the  data  in  Figures  1 5  and  1 6,  which 


assumed  a  linear  relationship  between  a  and  c,  yielded 
the  following  empirical  relationship: 


291  (1.55  -c)h 
h  +  2.95 


(21) 


where  c  is  expressed  in  percent,  h  in  cm  and  a  in  kPa. 


Urea  Concentration  {%) 

Figure  1 7.  Variation  of  flexural  strength  during 
freeze-up  with  urea  concentration  in  water  and 
ice  thickness. 
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Figure  18.  Effect  of  loading  direction  on  measured 
flexural  strength. 


Figure  19.  Variation  of  o' /o  with  relative  top  layer 
thickness  hf/h. 


The  effect  of  the  direction  of  load  application  on 
the  flexural  strength  calculated  from  cantilever  beam 
tests  according  to  eq  1  is  evidenced  in  Figure  18, 
where  o',  the  strength  obtained  with  upward  loading, 
i.e.  with  the  ice  bottom  in  tension,  is  plotted  against 
o,  the  strength  calculated  under  downward  loading, 
i.e.  with  the  top  surface  of  the  ice  in  tension.  In 
Figures  18  and  19  the  data  presented  were  obtained 
during  both  freeze-up  and  warm-up.  The  correspond¬ 
ing  values  of  o' jo  were  found  to  vary  from  as  low  as 
0.20  to  as  high  as  1 .0  (in  one  instance  only),  with 
80%  of  the  data  falling  between  1/3  and  2/3.  The 
values  of  o'/o  have  been  plotted  versus  ht/h,  the 
ratio  of  top  layer  thickness  of  total  ice  thickness,  in 
Figure  19,  where  lines  corresponding  to  eq  18  of  the 
proposed  two-layer  model  have  been  drawn  for  var¬ 
ious  values  of  f,  /E2  .  It  can  only  be  said  that  the 
experimental  data  exhibit  a  general  increase  of  o'/o 
with  hjh  as  predicted  by  the  model.  It  is  worth 


mentioning  that  measurements  performed  on  fresh¬ 
water  ice  also  showed  that  the  strength  o'  obtained 
by  upward  loading  is  also  about  80%  of  the  strength 
o  obtained  by  downward  loading  (see  e.g.  Timco  and 
Frederking  1981).  This  difference  is  usually  attrib¬ 
uted  to  a  difference  in  ice  crystal  size  between  the 
top  and  bottom  surfaces  of  the  ice.  This  is  character¬ 
istic  also  of  urea-doped  ice,  but  was  not  taken  into 
consideration  in  the  proposed  two-layer  model,  which 
assumes  uniform,  isotropic  properties  in  each  of  the 
two  layers. 

During  the  period  J uly-October  of  1981,  syste¬ 
matic  ice  strength  measurements  were  performed  at 
three  locations  along  the  south  wall  of  the  tank  and 
at  various  times  during  freeze-up  and  warm-up.  At 
each  of  the  three  locations,  four  cantilever  beam  tests 
were  carried  out  to  verify  the  repeatability  of  the 
results;  comparison  of  the  results  obtained  at  the 
three  locations  at  any  given  time  would  permit  a 
check  of  the  uniformity  of  ice  strength  in  the  tank. 

In  all,  137  sets  of  four  beam  tests  were  carried  out 
at  47  different  times.  The  results  of  these  tests  have 
been  listed  in  Appendixes  B  and  C.  For  each  set  of 
four  tests  at  one  location  and  at  one  time,  the  local 
mean  ice  strength,  local  standard  deviation  and  local 
coefficient  of  variation,  i.e.  the  ratio  of  standard 
deviation  to  mean  strength,  were  calculated.  Sim¬ 
ilarly,  for  each  of  the  47  test  times,  the  overall  mean 
strength,  overall  standard  deviation  and  overall  coef¬ 
ficient  of  variation  were  also  calculated.  The  results 
of  these  computations  are  listed  in  Table  4,  and  the 
histograms  of  the  local  and  overall  coefficients  of 
variation  are  shown  in  Figures  20  and  21 ,  respectively. 
The  mean  local  coefficient  of  variation  was  7.94% 
with  a  standard  deviation  of  5.65%,  and  the  overall 
coefficient  of  variation  was  6.31%  with  a  standard 
deviation  of  4.54%.  These  results  indicate  satisfactory 
repeatability  and  uniformity  in  the  ice  strength  meas¬ 
urements,  especially  when  compared  with  similar  re¬ 
sults  available  in  the  literature  on  ice  testing. 

Elastic  modulus 

The  elastic  modulus  of  an  ice  sheet  in  the  test  ba¬ 
sin  was  measured  by  the  plate  deflection  method.  In 
most  cases  the  load  was  applied  near  the  center  of 
the  basin.  In  one  instance,  the  uniformity  of  the  elas¬ 
tic  modulus  was  checked  by  repeating  the  measure¬ 
ments  at  18  different  locations  in  the  tank.  The  re¬ 
sults  of  this  check  are  given  in  Figure  22;  the  mean 
elastic  modulus  was  found  to  be  54.8  MPa,  with  a 
standard  deviation  of  5.2  MPa  or  about  10%  of  the 
mean,  which  was  considered  quite  satisfactory. 

The  effects  of  ice  thickness  and  urea  concentra¬ 
tion  on  the  elastic  modulus  of  urea-doped  ice  during 
freeze-up  are  presented  in  Figures  23  and  24.  Figure 
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Table  4.  Results  of  systematic  measurements  of  ice  bending  strength 


Date  and  time 

Local 

mean 

strength 

(kPa) 

Local 

strength 

deviation 

( kPa ) 

Local 
coef.  of 
variation 

(%> 

Overall 

mean 

strength 

(kPa) 

Overall 

strength 

deviation 

(kPa) 

Overall 
coef.  of 
variation 

(%) 

99.6 

6.5 

6.5 

8  July,  1300 

88.1 

4.0 

4.5 

90.7 

6.4 

7.1 

84.5 

3.8 

4.5 

83.8 

2.7 

3.2 

9  July,  0800 

81.1 

7.1 

8.8 

86.7 

6.2 

7.1 

95.3 

9.2 

9.6 

57.1 

2.1 

3.6 

9  July,  1600 

72.4 

4.3 

5.9 

69.7 

9.4 

13.5 

79.6 

3.7 

4.6 

181.0 

8.8 

4.9 

30  Aug,  0930 

174.0 

7.5 

4.3 

173.3 

6.6 

3.8 

165.0 

12.1 

7.3 

200.0 

11.9 

5.9 

30  Aug,  1 900 

200.0 

11.0 

5.5 

193.4 

9.3 

4.8 

180.3 

15.0 

8.3 

167.5 

13.7 

8.2 

31  Aug,  0845 

170.5 

8.9 

5.2 

167.1 

3.0 

1.8 

163.2 

1.8 

1.1 

148.2 

10.8 

7.3 

31  Aug,  1700 

143.3 

2.5 

1.7 

145.3 

2.1 

1.4 

144.5 

11.2 

7.7 

99.5 

2.9 

3.0 

31  Aug,  2230 

101.6 

4.3 

4.2 

100.6 

0.85 

0.8 

100^6 

1.6 

1.6 

75.2 

2.8 

3.8 

1  Sept,  0830 

69.8 

5.0 

7.2 

71.3 

2.8 

3.9 

69.0 

3.3 

4.8 

60.9 

1.6 

2.6 

1  Sept,  1530 

60.4 

2.6 

4.3 

59.0 

2.3 

4.0 

55.7 

6.1 

10.9 

51.9 

2.0 

3.8 

1  Sept,  2230 

52.6 

0.7 

1.2 

5  3.2 

3.6 

6.7 

48.1 

1.0 

2.0 

2  Sept,  0840 

35.8 

19.4 

54.2 

44.5 

6.2 

13.9 

49.6 

S.l 

10.3 

140.9 

11.4 

8.1 

6  Sept,  1015 

132.6 

8.6 

6.5 

137.7 

3.7 

2.7 

139.7 

5.1 

3.6 

148.6 

4.0 

2.7 

6  Sept,  1800 

124.7 

12.3 

9.9 

133.8 

10.6 

7.9 

128.0 

5.1 

4.0 

106.8 

14.7 

13.8 

7  Sept,  1630 

101.0 

10.0 

9.9 

101.0 

4.7 

4.7 

95.2 

7.9 

8.3 

75.6 

1.7 

2.3 

7  Sept,  2200 

80.3 

1.6 

2.0 

77.2 

2.2 

2.8 

75.8 

1.4 

1.8 

44.7 

3.6 

8.0 

8  Sept,  0830 

40.5 

6.0 

14.8 

43.3 

2.0 

4.6 

44.8 

5.9 

13.1 
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Local 

Local 

Local 

Overall 

Overall 

Overall 

mean 

strength 

coef.  of 

mean 

strength 

coef.  of 

strength 

deviation 

variation 

strength 

deviation 

variation 

Date  and  time 

(kPa) 

(kPa) 

(kPa) 

(kPa) 

(%) 

118.5 

2.6 

2.2 

10  Sept,  2130 

85.0 

11.9 

14.0 

102.2 

13.7 

13.4 

103.0 

10.0 

9.7 

126.1 

9.0 

7.1 

1 1  Sept,  0900 

118.2 

9.5 

8.0 

125.0 

5.2 

4.1 

130.7 

3.4 

2.6 

105.1 

3.0 

2.9 

11  Sept,  1230 

91.1 

4.6 

5.3 

98.9 

5.8 

5.9 

100.6 

5.7 

5.6 

75.6 

3.5 

4.6 

11  Sept,  1700 

67.4 

9.7 

14.3 

71.8 

3.4 

4.7 

72.5 

2.5 

3.5 

40.4 

4.3 

10.6 

12  Sept,  0930 

29.0 

2.1 

7.r 

32.0 

6.0 

18.7 

26.7 

2.8 

10.4 

28.3 

3.0 

10.4 

12  Sept,  1945 

30.8 

1.6 

5.2 

28.8 

1.5 

5.1 

27.3 

1.5 

5.4 

81.3 

8.7 

10.6 

15  Sept,  0930 

88.8 

2.8 

3.2 

87.1 

4.2 

4.8 

91.2 

9.8 

10.7 

59.6 

14.7 

24.7 

15  Sept,  1145 

62.8 

2.4 

3.8 

62.3 

2.0 

3.2 

64.4 

2.4 

3.7 

37.3 

2.6 

7.0 

15  Sept,  1715 

38.0 

7.5 

19.8 

38.1 

0.74 

1.9 

39.1 

5.9 

15.0 

23.4 

4.9 

20.7 

15  Sept,  2400 

20.4 

1.6 

8.0 

22.8 

1.7 

7.6 

24.5 

3.3 

13.4 

17.9 

0.7 

3.8 

16  Sept,  0930 

14.2 

1.6 

11.2 

15.6 

1.6 

10.4 

14.8 

2.1 

14.1 

61.8 

3.6 

5.8 

20  Sept,  0815 

69.6 

3.6 

5.2 

65.1 

3.3 

5.1 

63.9 

3.9 

6.1 

46.2 

4.3 

9.2 

20  Sept,  1430 

42.5 

4.8 

11.3 

42.5 

3.0 

7.1 

38.8 

3.2 

8.3 

29.0 

4.0 

13.7 

20  Sept,  2200 

28.6 

5.2 

18.0 

30.5 

2.5 

8.0 

34.0 

4.6 

13.4 

13.5 

2.9 

21.5 

21  Sept,  0930 

17.7 

1.3 

7.4 

13.7 

3.2 

23.3 

9,9 

0.9 

8.9 

67.5 

3.6 

5.3 

27  Sept,  0900 

62.4 

2.0 

3.3 

66. 5 

3.1 

4.6 

69.7 

5.9 

8.5 

27  Sept,  1345 

39.7 

43.5 

11.4 

5.5 

28.7 

12.7 

41.6 

1.9 

4.6 

27  Sept,  1815 

21.8 

23.3 

2.7 

2.3 

12.2 

9.7 

22.6 

0.8 

3.3 

4 
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Frequency 


Date  and  firm 

Local 

mean 

strength 

fkPa) 

Local 

strength 

deviation 

(kPa) 

Local 
coef.  of 
variation 
(kPa) 

Overall 

mean 

strength 

(kPa) 

Overall 

strength 

deviation 

(kPa) 

O  verall 
coef.  of 
variation 
(kPaj 

80.6 

4.7 

5.9 

7  Oct,  1300 

74.3 

3.7 

5.0 

76.4 

3.0 

3.9 

74.3 

1.2 

1.6 

127.0 

7.4 

5.8 

8  Oct,  0850 

124.6 

12.0 

9.6 

122.5 

4.8 

3.9 

115.8 

6.6 

6.6 

10.4 

2.2 

2.4 

8  Oct,  1230 

94.7 

1.6 

9.6 

89.9 

4.2 

4.6 

84.5 

6.1 

6.6 

S8.8 

3.6 

6.1 

8  Oct,  1730 

60.1 

0.5 

0.8 

61.5 

3.0 

4.8 

65.6 

4.2 

6.5 

105.2 

14.5 

13.8 

9  Oct,  0945 

98.4 

4.7 

4.8 

101.0 

3.0 

3.0 

99.3 

6.9 

6.9 

66.5 

3.1 

4.7 

10  Oct,  0915 

76.8 

2.5 

3.2 

70.8 

4.4 

6.2 

69.1 

1.8 

2.6 

50.8 

6.3 

12.3 

11  Oct,  1145 

46.7 

4.15 

8.9 

42.5 

2.1 

4.9 

86.9 

7.0 

8.1 

17  Oct,  2120 

91.5 

21.1 

23.6 

81.6 

10.9 

13.3 

66.5 

3.2 

4.9 

115.0 

8.3 

7.3 

18  Oct,  1410 

99.2 

3.6 

3.6 

106.9 

6.4 

6.0 

106.6 

4.4 

4.2 

104.6 

8.2 

7.8 

19  Oct,  1650 

96.3 

6.9 

7.2 

96.5 

6.6 

6.8 

88.5 

3.2 

3.6 

63.6 

7.9 

12.4 

20  Oct,  0845 

52.0 

8.1 

15.5 

55.4 

5.8 

10.5 

50.7 

5.0 

9.8 

49.2 

3.9 

7.9 

20  Oct,  1545 

46.1 

2.5 

5.4 

47.2 

1.4 

3.0 

46.3 

8.8 

19.1 

v, Coefficient  of  Voriotion  (%) 


v,  Coefficient  of  Variation  (%) 


Figure  20.  Local  coefficient  of  variation  of  flex-  Figure  21 .  Overall  coefficient  of  variation  of  flex¬ 
ural  strength.  ural  strength. 


(m)  34  32  30  28  26  24  22  20  1 8  16  14  1 2  10  0  6  4  2  0 


Figure  22.  Distribution  of  elas tic  modulus  in  the  tank  (22  Nov  1 981  j. 


0  12  3  4  5  6 


Ico  Thickness  (cm) 

Figure  23.  Effect  of  ice  thickness  on  elastic  modulus. 


23  shows  that,  except  for  the  results  for  thin  ice 
sheets  (of  thickness  of  2  to  2.5  cm)  where  the  loads 
which  can  be  applied  are  very  small  and  accuracy  in 
the  measurements  is  questionable,  the  elastic  modulus 
is  practically  independent  of  ice  thickness.  On  the 
other  hand,  the  elastic  modulus  decreases  rapidly 
with  increasing  urea  concentration  from  about  300 
MPa  at  c  =  0.45%  to  approximately  70  MPa  at  c  = 
1.05%. 


0  SO  100  ISO  200  250 

E,  Downward  Elastic  Modulus  (MPo) 


Figure  25.  Effect  of  loading  direction  on 
elastic  modulus. 
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Concentration  o<  Urea  Solution  (%) 


Figure  24.  Effect  of  urea  concentration  on 
elastic  modulus. 


Finally,  as  shown  in  Figure  25,  which  includes 
measurements  taken  during  freeze-up  and  warm-up 
of  the  ice  sheets,  the  elastic  modulus  was  found  to  be 
independent  of  the  direction  (upward  or  downward) 
of  load  application. 

Because  of  insufficient  data,  the  effect  of  temper¬ 
ature,  if  any,  on  the  elastic  modulus  could  not  be  in¬ 
vestigated  with  any  degree  of  confidence. 

Ratio  E/o 

The  value  of  E,  elastic  modulus,  is  plotted  against 
a,  the  flexural  strength  in  the  downward  mode  of 
load  application,  in  Figure  26.  All  the  data  are  seen 
to  lie  between  the  two  lines  El  a  =  500  and  E/o  = 
2000.  The  ratio  E/o  itself  was  plotted  versus  urea 
concentration  c  in  Figure  27.  Figure  27  indicates 
that  E/o  decreases  with  urea  concentration  but 
appears  to  approach  a  constant  value  of  about  1000 
for  c  equal  to  0.90%  or  more.  The  average  values  of 
E/o  are 

E/o  =  1 800  ±  300  for  c  =  0.45% 

E/o=  1250  ±  350  fore  =  0.70% 

E/o=  1000  ±  250  fore  =  0.95% 

E/o  =  970  ±300  for  c  =  1 .05% 


The  decrease  of  E/a  with  increasing  solute  con¬ 
centration  in  the  meltwater  was  also  observed  for 
saline  ice  (Schwarz  1975),  and  predicted  analytically 
by  Weeks  and  Assur  (1967),  who  derived  a  relation¬ 
ship  between  E/a  and  brine  volume  in  the  ice. 

Finally,  in  Figure  28  the  ratio  E/a  is  plotted 
against  hjh.  In  spite  of  the  scatter  in  the  data  points, 
it  can  be  seen  that  E/a  generally  increases  with  in¬ 
creasing  hjh,  at  least  in  the  range  of  this  latter  pa¬ 
rameter  for  the  present  study.  This  variation  of  E/a 
with  hjh  is  in  qualitative  agreement  with  the  two- 
layer  model  prediction  as  illustrated  by  the  curve 
drawn  in  Figure  28  which  corresponds  to  eq  19  with 
arbitrarily  selected  values  of  EJol  =  4000  and  F,  /F3 
=  100.  Since  no  ice  sheet  with  values  of  hjh  less 
than  0.1  was  grown  in  the  tank,  verification  of  the 
reverse  behavior  of  E/a  with  hjh  for  small  hjh  was 
not  possible.  However,  Timco  (pers.  comm.)  has 


Figure  28.  Variation  of  E/o  with  rel¬ 
ative  top  layer  thickness  h/h. 


Figure  2  7.  Variation  of  E/o  with  urea 
concentration. 


indicated  that  the  urea-doped  ice  sheets  grown  in  the 
test  tank  of  the  National  Research  Council  of  Canada 
exhibit  a  very  thin  top  layer,  and  that  values  of  E/a 
of  1 500  to  2000  were  commonly  measured.  It  is 
likely  that  the  difference  in  top  layer  thickness  be¬ 
tween  the  two  tanks  is  due  to  a  difference  in  growth 
rate  (higher  in  the  NRC  tank  where  cold  air  is  blown 
over  the  tank). 

Properties  of  urea-doped  ice  during  warm-up 

As  mentioned  in  the  Introduction ,  tempering  of 
the  ice  sheet  is  the  practical  method  for  achieving 
the  flexural  strength  values  required  in  most  mode) 
tests,  which  are  usually  lower  than  those  obtained 
during  freeze-up. 

As  soon  as  the  refrigeration  system  in  the  tank 
room  is  shut  off,  the  air  temperature  rises  naturally, 
as  shown  in  Figure  29,  with  corresponding  variation 
in  the  ice  properties,  as  follows. 

Flexural  strength 

As  the  air  temperature  in  the  room  rises,  so  does 
the  ice  surface  temperature;  consequently  the  strength 
of  the  sheet’s  top  layer  and,  therefore,  the  nominal 
strength  of  the  ice  sheet,  decreases.  Additional  de¬ 
creases  in  the  ice  strength  can  be  related  to  enlarge¬ 
ment  of  the  brine  pockets  in  both  the  top  and  bottom 
layers.  As  previously  mentioned,  Weeks  and  Assur 
(1967)  showed  that  the  strength  of  sea  ice  decreases 
with  brine  volume,  which  itself  increases  with  sur¬ 
rounding  air  temperature. 

Examples  of  strength  variation  with  time  are 
shown  in  Figure  30  for  three  ice  sheets  grown  from 
the  same  0.95%  urea  solution  and  with  nearly  identi¬ 
cal  thicknesses  and  initial  strength  at  the  start  of  the 
warm-up  period.  The  final  strength  for  the  tests  was 
targeted  to  20  kPa. 
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Temperoture 

(°C)|  (°F)[ 


2  4  ("fO 

Warm-up  T ime 

Figure  29.  Example  of  air  temperature  change  during  warm-up. 


I  j-  Time  when 
r  warm  -up  Starts 


In  Figure  31 ,  the  data  collected  for  the  ice  sheets 
grown  for  the  four  urea  concentrations  investigated 
are  presented  as  o/oi  versus  t/hf,  where  o  is  the  nom 
inal  ice  strength  at  time  t  during  warm-up,  o;  is  the 
initial  strength  at  t  =  0,  and  ft,  is  the  final  desired  ice 
thickness.  The  parameter  t/h? ,  where  hi  is  the  ice 
thickness  at  the  start  of  warm-up,  t  =  0,  could  also 
have  been  used  since  it  was  noted  that,  for  the  range 
of  experimental  conditions  studied  here,  h{  and  h{ 
are  proportional  to  each  other  (see  Fig.  8). 


Figure  30.  Examples  of  decrease  in  flexural  strength 
during  warm-up. 


•  May  6  -  8  (3.2) 

a  Jul.  7  -  II  (5.4) 

o  Jul  22-25  (4.2! 


h,(cm) 

h((cm 

A  Sep 26-28 

(3.1) 

o  Nov  19-23  (No  31(5  6) 

A  Oct  6-12 

(5.  1) 

V  Nov  25-27  (No4)(3 .4) 

o  Oct  16-21 

(5.0) 

<j  Nov27-Dec  l(No.5H5.6) 

•  Nov  12-16  (No.  1  )  (5  6) 

T  Dec  1-3  (N0.6X34) 

■  Nov  17-19  (N0.2H33) 

A  \  v 

06  —  \,  V 


1.0  1.5 


a.  c  =  1.05% 


Figure  31.  Plot  of  o/o,  vs  t/h?. 


b.  c  =  0.95% 


Figure  31  (cont d). 


Figure  32.  Plot  of  o/ot  vs  t/h7  for  all  urea  concentrations. 


In  Figure  31 ,  which  shows  that  the  data  for  each 
urea  concentration  collapse  reasonably  well,  empiri¬ 
cal  relations  of  the  form 

£-  =  i  -p(t/hff 

ai 

have  been  determined  by  nonlinear  regression  anal¬ 
ysis,  namely 


c  =  1 .05% 

a/Oj  =  1  -0.70  (tlh7fA* 

(22a) 

c  =  0.95% 

o/Oj  =  1  -0.64  (tfhj)oso 

(22b) 

c  =  0.70% 

o/Oj  =  1  -0 .67  (f/Zif)0*33 

(22c) 

c  =  0.45% 

a/Oj  =  1  -  0.70  (t/h7)033 

(22d) 

The  data  for  the  f  our  urea  concentrations  have 
been  collected  on  the  single  plot  shown  in  Figure  32, 
from  which  the  following  empirical  relations  have 
been  calculated: 

jj-  =  1  -  (0.25  t/h7')0-33  (23) 

and 

~=  exp  [-(1.1 2  t/h7)063}  (24) 

In  eq  24  the  factor  K  =  1.12  cm2/hr  might  be  related 
to  the  thermal  diffusivity  of  the  urea-doped  ice.  From 
Figure  32,  it  can  be  seen  that  eq  23  is  a  good  approx¬ 
imation  to  eq  24  for  the  range  0  <  t/h7  <  3. 
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O  I  2 

Vh*  (hr/cm2) 

Figure  33.  Plot  of  E/E.  vs  t/h?  /or  all  urea  concentrations. 


Elastic  modulus 

The  ratio  £/£,,  where  Ei  is  the  elastic  modulus  at 
the  start  of  warm-up,  has  been  plotted  against  t/h* 
in  Figure  33  for  all  four  urea  concentrations. 

Two  empirical  relationships  have  been  determined 
by  nonlinear  regression  analysis: 


§7  =  1 -(0.28  t/hj  )° 


=exp  [-(1.63  t/hf)0"54]  .  (26) 

ti 

These  two  relationships  have  also  been  plotted  in 
Figure  33. 

Ratio  E/o 

Comparison  of  Figures  32  and  33  indicates  that 
the  elastic  modulus  E  decays  slightly  faster  than  the 


flexural  strength  o  with  time  during  warm-up.  There 
fore,  some  decrease  in  the  ratio  E/o  may  be  expected 
during  warm-up.  However,  because  of  the  experi¬ 
mental  scatter  in  the  data  for  both  o  and  E,  even 
larger  scatter  may  be  expected  in  their  ratio.  The 
quantity  E*  =  (E/o)/(EJof)  has  been  plotted  versus 
t/hj  in  Figure  34.  It  can  be  seen  that  40  out  of  49 
points,  or  80%  of  the  data,  falll  between  E*  =  0.55 
and  1 .20,  and  that  the  rather  considerable  scatter  in 
the  data  points  does  not  allow  identification  of  any 
significant  variation  of  E*  with  t/hj . 

The  measurements  of  E  and  a  are  also  collected 
in  Figure  35  as  E  versus  a  for  the  four  urea  concen¬ 
trations  used  in  the  study.  From  these  figures  it 
appears  that  the  majority  of  the  data  fall  between 
the  lines  E/o  =  500  and  E/o  =  1 500,  and  that  the 
ratio  E/o  generally  decreases  with  increasing  urea 
concentration,  as  mentioned  previously. 


E,  Elcutic  Modulot  (MPo)  E-  Eloi"c  Woduiui  U*Po> 


Figure  35.  Plot  of  E  vs  o. 


APPLICATIONS  TO  TEST  PROGRAM 
SCHEDULING 

Any  test  of  mechanical  structures  in  level  ice  can 
be  conducted  only  when  the  required  ice  character¬ 
istics,  namely  thickness  and  flexural  strength,  are 
achieved  in  the  tank;  in  other  words  the  ice  controls 
the  time  at  which  model  tests  are  performed.  The 
results  of  the  present  study  have  been  used  to  form¬ 
ulate  a  reliable  time  schedule  for  growing  and  temper¬ 
ing  an  ice  sheet  to  reach  present  conditions,  so  that 
tests  can  be  started  in  the  morning  at  regular  working 
hours  rather  than  at  odd  hours  of  the  day  or  night. 

The  method  detailed  below  is  valid  for  the  0.95% 
urea  concentration  now  being  used  in  the  CRREL 
test  basin,  and  for  the  particular  refrigeration  cap¬ 
abilities  of  the  CRREL  plant. 


From  eq  20b,  22b  and  6  a  family  of  curves  have 
been  plotted  in  Figure  36  for  a  wide  range  of  ice 
thicknesses,  flexural  strengths,  and  warm-up  times. 
When  the  flexural  strength  and  ice  thickness  are  given 
as  test  conditions  to  be  achieved  in  the  tank,  the  cor¬ 
responding  initial  ice  thickness  and  necessary  warm¬ 
up  time  can  be  determined  from  Figure  36.  Once  the 
initial  ice  thickness  is  determined,  the  necessary  freeze- 
up  time  and  temperature  are  obtained  from  Figure 
37,  where  eq  4  has  been  plotted  as  h{  versus  time  for 
several  air  temperatures,  with  the  values  of  the  coef¬ 
ficients  k{  and  hja  obtained  from  Figure  6b.  For 
example,  to  obtain  a  model  ice  sheet  of  final  thick¬ 
ness  hf  =  4.5  cm  and  final  flexural  strength  o  =  50 
kPa,  Figure  36  indicates  that  warm-up  should  be 
started  when  the  ice  thickness  reaches  a  value  hi  = 

3.9  cm  and  should  last  for  approximately  18  hours. 
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Figure  36.  Diagram  for  planning  model  ice  prep¬ 
aration  (c  =  0.95%). 


Then  from  Figure  37,  the  freeze-up  time  required 
to  grow  a  3.9-cm-thick  ice  sheet  is  determined  as  22 
hours  at  an  air  temperature  of  -1 8°C  (0°F)  or  32 
hours  at  an  air  temperature  of  -1 2°C  (10°F). 

Once  the  total  time  (freeze-up  and  warm-up)  to 
achieve  required  ice  conditions  has  been  determined, 
either  the  time  at  which  seeding  should  take  place 
can  be  decided,  or  the  time  at  which  the  test  is  likely 
to  be  performed  can  be  calculated  with  a  good  degree 
of  confidence.  In  the  above  example,  the  total  time 
(assuming  a  -12°C  air  temperature  during  freeze-up) 
is  50  hours.  Therefore,  if  the  test  is  scheduled  for  a 
Monday  morning  at  0800  hours,  the  tank  should  be 
seeded  the  preceding  Saturday  morning  at  about  0600 
hours.  On  the  other  hand,  if  for  some  reason  seeding 
cannot  be  done  before  1 200  hours  on  Saturday,  the 
start  of  the  test  can  be  scheduled  for  Monday  after¬ 
noon  at  about  1400  hours,  unless  the  air  temperature 
in  the  room  during  freeze-up  can  be  adjusted  to  -15°C, 
in  which  case  the  test  is  likely  to  be  possible  in  the 
early  morning. 

This  scheduling  method  has  been  in  operation 
now  for  some  time  at  CRREL  and  has  proven  to  be 
reliable.  For  example,  the  tests  whose  results  are 
shown  in  Figure  30  were  so  scheduled. 

It  goes  without  saying  that  although  this  method 
can  and  does  reduce  the  time  uncertainty  in  sched¬ 
uling  tests  in  ice,  monitoring  of  air  temperature  and 
ice  characteristics  (thickness  and  strength)  remains 
necessary,  but  usually  less  frequently,  to  account 
for  any  variation  in  plant  operation  or  outside  air 
temperature,  and  to  allow  for  unforeseen  circum¬ 
stances. 


Figure  37.  Expected  ice  thickness  vs  freeze-up 
time  for  different  air  temperatures. 

SUMMARY  AND  CONCLUSIONS 

A  study  of  the  growth  process  and  mechanical 
properties  of  some  50  urea-doped  ice  sheets  in  the 
CRREL  ice  test  basin  has  led  to  the  following  results: 

1 .  Uniformity  of  ice  thickness  in  the  main  test 
area  of  the  basin  is  considered  satisfactory,  with  a 
standard  deviation  in  the  ice  thickness  of  1  mm  or 
less  for  thicknesses  of  2  to  7  cm. 

2.  A  relationship  between  ice  thickness  and  neg¬ 
ative  degree-hours  was  well  established. 

3.  The  structure  of  urea-doped  ice  was  found  to 
be  similar  to  that  of  sea  ice  except  for  a  relatively 
thick  top  layer  over  a  bottom  layer  of  columnar 
structure  typical  of  dendritic  growth. 

4.  A  mathematical  two-layer  ice  model  was  de¬ 
veloped  to  predict  the  behavior  of  ice  properties  for 
comparison  with  the  results  of  laboratory  measure¬ 
ments.  There  is  qualitative  agreement  between  theory 
and  experiments.  Before  the  mathematical  two-layer 
model  can  be  improved,  better  understanding  of  the 
growth  mechanism  of  the  top  ice  layer  is  necessary. 

5.  Uniformity  of  ice  flexural  strength  and  elastic 
modulus  in  the  tank  was  found  to  be  satisfactory. 

6.  Empirical  formulas  were  obtained  that  give  the 
relationships  between  1 )  ice  thickness,  flexural  strength 
and  elastic  modulus  and  2)  air  temperature,  urea  con¬ 
centration,  and  ratio  of  top  layer  thickness  to  total 

ice  thickness  during  the  freeze-up  period. 

7.  During  the  tempering  or  warm-up  period,  the 
reduction  in  strength  and  in  elastic  modulus,  and  the 
increase  in  ice  thickness,  were  quantified,  and  em¬ 
pirical  relationships  were  established. 
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8.  The  ratio  E/a  was  found  to  vary  widely  during 
both  freeze-up  and  warm-up  from  one  ice  sheet  to 
another  grown  under  similar  conditions.  On  the 
average,  this  ratio  could  be  larger  than  2000  for  ice 
sheets  grown  from  a  0.45%  urea  solution  in  water, 
and  smaller  than  1000  for  ice  sheets  grown  from  a 

1 .05%  solution.  The  two-layer  model  prediction  that 
a  considerable  increase  in  E/a  could  be  achieved  with 
ice  sheets  with  a  thin  top  layer  (less  than  5%  of  the 
total  thickness)  could  not  be  verified. 

9.  The  results  of  the  study  were  combined  into 
charts  that  permit  establishment  of  a  reliable  test 
schedule  for  producing,  in  the  tank,  ice  sheets  with 
the  required  thickness  and  flexural  strength.  These 
charts  have  been  used  successfully. 

In  addition  to  the  above  results,  the  study  has 
confirmed  that  urea-doped  ice  is  a  very  satisfactory 
alternative  to  saline  ice  as  model  ice  for  testing  ice- 
structure  interaction.  Urea-doped  ice  allows  a  model 
scale  as  low  as  1  /40,  and  practically  eliminates  cor¬ 
rosion  problems  in  experimental  facility  equipment. 

rhe  study  has  also  led  to  improvements  in  the 
operation  of  the  ice  tank  at  CRREL,  and  has  indi¬ 
cated  what  further  modifications  and  improvements 
are  necessary  or  would  be  useful.  In  particular,  meas¬ 
urement  techniques  for  both  flexural  strength  and 
elastic  modulus  could  be  improved,  and  application 
of  different  techniques  would  be  useful  for  com¬ 
parison  between  laboratory  and  field  data.  Further 
research  efforts  are  needed  to  clarify  the  growth 
process  of  urea  ice  and  its  two-layer  structure,  in 
order  to  be  able  to  control  the  thickness  of  the  top 
layer  and  ach.  ve  as  high  values  of  the  ratio  E/a  as 
possible. 
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APPENDIX  A:  RESULTS  OF  ICE  THICKNESS  MEASUREMENTS 
FOR  VARIOUS  GROWTH  CONDITIONS 
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-435.1 

** 

42.5 

0 

OO 

1  ,C8 

-532.5 

1  -  2  Dec 

(seeded  1700  1  Dec) 


15.5 

20.0 


2.21 

2.76 


0.6  5 
0.73 


■222.6 

•293.5 


Note:  vrea  concentration  c  =  0.70?  f-om  5  to  21  Sept 

c  =  0.95?  from  2-  Sept  on 


i 


i 
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APPENDIX  B:  PROPERTIES  OF  UNTEMPERED  ICE 


Experiment 

'Jrea 

T 1  -ne 

3  earn 

F 1 exura 1 

F.  1  as  1 1  c 

p  er 1 od 

concentrat 1  on 

of  beam  tests 

th  lekness 

Strength 

modu 1 u  s 

(%) 

(cm) 

(kPa) 

(MPa) 

27-29  May  1981 

1.05 

28  May  - 

0930 

2.5 

42.2 

t  6.1 

2.6 

45.3 

i  2.0 

28  May  - 

1500 

3.1 

34.3 

t  0.2 

2.9 

35.1 

t  1.3 

23  May  - 

2100 

3.3 

24.3 

t  2.3 

3.1 

39.3 

t  4.9 

29  May  - 

0800 

3.4 

24.1 

t  2.1 

3.5 

21.9 

t  1.4 

29  May  - 

1400 

3.65 

20.3 

t  1.3 

31.  < 

1-3  June  1981 

1.05 

2  June  - 

0900 

1.75 

56.4 

t  6.4 

1.9 

58.6 

t  5.9 

1.9 

37.9 

t  16.8 

2  June  - 

1600 

2.4 

93.2 

t  12.0 

2.4 

92.0 

t  4.6 

2.4 

66.6 

t  4.0 

3  June  - 

0900 

5.75 

72.3 

t  4.5 

4.1 

79.2 

t  5.0 

3.9 

100 

i  12.3 

3  June  - 

1030 

3.9 

1 10 

t  3.2 

4.0 

75 

±  4.1 

3  June  - 

1200 

4.0 

87.7 

t  5.3 

4.2 

84.7 

t  3.2 

4.1 

32.8 

±  3.8 

25!.' 

3-5  June  1981 

1.05 

4  June  - 

1030 

1.85 

68.4 

±  4.4 

2.0 

60.3 

±  3.9 

91. 

1.9 

71.2 

±  3.1 

4  June  - 

16  30 

2.3 

80.9 

±  1.4 

2.45 

70.9 

(M 

• 

+1 

2.4 

82.9 

±  8.1 

5  June  - 

0830 

3.8 

78.2 

±  2.3 

3.9 

89.0 

±  5.7 

3.75 

94.1 

±  10.2 

5  June  - 

1  100 

4.0 

103.4 

+  11.2 

4.1 

87.1 

±  3.3 

23.2 

4.1 

11 1 .4 

+  9^ 

3! 


Experiment  Urea  Time  Beam  Flexural  Elastic 

period  concentration  of  beam  tests  Ihlckness  strenqlh  modulus 

<*)  (cm)  (kPa)  (MPa) 


5-8  June  >981 


8-11  June  1981 


6  June 

6  June 

7  June 

7  June 

8  June 
8  June 


10  June 


1.7 

47.5  ±  15.9 

1.8 

52.9  ±  2.9 

1.9 

48.8  ±  3.2 

2.85 

59.4  +  3.8 

2.9 

69.0  ±  13.9 

5.1 

62.3  ±  1.9 

3.75 

70.7  +  3  Jc> 

3.85 

70.5  +  10.4 

4.1 

D8.8  ±  4.5 

4.4 

75.0  ±  5.5 

4.6 

72 j6  ±  3.6 

4.8 

53.3  ±  3.2 

5.2 

94.0  +  4.5 

5.6 

85.5  +6.1 

5.4 

77.3  ±  2.2 

5.9 

6  2.1  +  7.8 

10  June 


11  June 


81.6  ±  7.2 

74.7  t  7.5 

57.7  ±  5.6 

83.1  +  8.1 

76.5  ±  4.9 
79.4  +  7.1 

67.2  ±  3.6 

91.9  +  3.8 

89.8  ±4.6 

90.8  +6.9 
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Urea  Time  Beam  Flexural  Elastic 

concentration  of  beam  tests  Ihlckness  strength  modulus 
(<) _ (cm) _ (kPa  ) _ (MPa) 

3.45  58.3  *  1.4 

3.4  58.2  ±  3.2 

99.6  *  8.3 
88.9  *  14.3 
90.0  *6.2 

95.8  *  7.9 _ 

19-22  June  1981  1.05  20  June  -  0800  2.1  67.7  ±  5.8 

2.2  66.4  ±  7.0 

2.4  70.4  ±4.6 

21  June  -  0800  4.4  87.6  *  3.8 

4.5  91.4  ±  8.6 

4.7  90.3  ±6.0 


22  June  -  0900  6.45  92.3  +  8.9 

_ 656 _ 89.2  ±  8.8 


24-25  June 

1981 

1.05 

25  June 

-  0830 

2.2 

55.5 

±  8.4 

2.4 

65.2 

i  3.9 

2.3 

68.1 

±  4.6 

26  June 

-  0800 

4.15 

95.3 

±  4.1 

4.25 

91.8 

±  4.2 

4.6 

88  JS 

±  5.4 

26-29  June 

1981 

1.05 

27  June 

-  0800 

1.9 

55.0 

±6.6 

2.1 

50.3 

±6.2 

2.0 

51.5 

±  1.1 

28  June 

-  0800 

4.0 

79.3 

±  7.0 

4.3 

81.0 

±  5.8 

4.1 

82.8 

±  3.5 

29  June 

-  0800 

6.25 

99.5 

±  7.0 

5.9 

77  JS 

±  18.3 

6.35 

92.0 

±  4.9 

29  June 

-  11  VI 

6.5 

73.3 

±  1.3 

6.2 

63.4 

±  0.9 

6.35 

93.2 

+  5.2 

30  June  -  1 

July 

1.05 

1  July  -  0800 

3.55 

83.5 

*  3.2 

1981 

3.7 

87.9 

±  7.2 

3.7 

88,6 

±6.0 

2-0  July  1981  1.05  3  July  -  1130  2.3  52.6  ±4.0 

2.5  63.0  ±  5.0 

2.4  58.5  ±  9.2 


Experiment 

period 
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i  -  * 

>v. 

Experiment 

Urea 

Time 

Beam 

F I exura 1 

E  lastic 

r-.- 

period 

concentration 

of  beam  tes  ts 

1h  Ickness 

strenqlh 

modulus 

(?) 

(cm ) 

(kPa) 

(MPa) 

70.6  +  2.2 

79.5  t  4.6 
74.3  +  5.3 

94.8  +  5.1 

91.2  +  27,8 

96.9  ±  12.2 

87.7  ±  3.2 

105.2  +6.8 
139.8  t  12.4 

124.5  +  18.9 

82.2  +  5.7 

106.3  +  9.2 
96.1  +  7.5 

83.6  +  1.6 

78.6  .+  9.2 

78.5  +  10.2 


13-15  July  1981  1.05 


15  July  -  0800 


95.9  +  8.0 
83.6  +  5.5 
102.5  t  7.2 


15  July  -  0930 


101.8  +  4.5 
86.4  +  9.3 
107.6  +  9.2 


19-21  July  1981  1.05 


20  July  -  0900 


87.0  +  2.6 
79.9  ±  5.2 
85.8  +  2  j6 


21  July  -  0800 


92.0  +  2.8 
87.5  +  4.1 
85.1  ±  4.7 


APPENDIX  C:  PROPERTIES  OF  TEMPERED  ICE 


Exp.  Per  lod  6-8  May  '81  Final  Ice  1tilckness  hf  *  3.2  cm 

Concentration  1.05)1 

Note:  Warm-up  started  7  May,  at  1230 

Flexural  Elastic  Warm- 


o 

Oate,  Ice  Itilckness  strenqlti  (kPa)  a/a 1  modulus  (MPa)  up  time  t/h 

time  h(cm)  a  a'  O'/a  <0|  *  75.2  kPa)  E  E'  t(hr)  f 


7  May 
1100 

2.63  -  2.78 

75.2  +  2.6 

7  May 

2.96 

58.5  ±  1.0 

0.78 

1.5 

0.147 

1400 

7  May 

3.20 

38.1  ±  4.1 

0.51 

3.5 

0.342 

1610 

7  May 

3.25 

20.8  ±  2.2 

0.21 

9.0 

0.879 

2130 

8  May 

3.20 

12.6 

0.17 

20.0 

1.953 

0815 

8  May 

3.10 

17.3  ±  1.0 

0.23 

21.5 

2.100 

1000 

8  May 

3.10 

15.0  ±  1.5 

0.20 

23.5 

2.294 

1200 

Exp.  Period  7-11  July 

•81 

Final  Ice  Itilckness  hf  * 

5.4  cm 

Concentration  1  ,05< 

Note: 

Warm-up  started 

9  July,  at  1030 

Flexural 

Elastic 

Warm- 

Date, 

Ice  Iti  Ickness 

strenqtt  (kPa) 

a/a  | 

modulus  (MPa) 

up  time 

t/ti  2 

time 

h  (cm) 

o  a' 

a1 /a  ((J|  *  84.7  kPa) 

E  E' 

t(hr ) 

f 

8  July 

2.9  -  2,95 

99.6  ±  6  .5 

1300 

3.0 

88.1  +  4.0 

60 

3.0 

84.5  ±  3.8 

9  July 

4.8 

83.8  ±  2.7 

0800 

4.8 

81.1  ±  7.1 

94 

4J55  -  4.7 

95.3  ±  9.1 

9  July 

5.1 

83.1  ±6.1 

0.89 

1300 

5.2 

89.8  ±  3.5 

0.96 

53 

2.5 

0.09 

5.2 

95.1  ±  7.7 

1.02 

9  July 

5.25 

57.1  ±  2.1 

0.6  1 

1600 

5.4 

72.4  ±  4.3 

0.78 

48 

6 

0.21 

5.3 

79J5  ±  3.7 

0.85 

35 


Date, 

time 

Ice  1h  ickness 
h  (cm) 

Flexural 
strenglh  (kPa) 

a  o' 

a'/a  <0| 

Elastic 

a/ at  modulus  (MPa) 

-  Ill  kPa )  E  E' 

Warm¬ 
up  time 
t(hr ) 

t/h2 

f 

9  July 

5.45 

64.8 

±  2.3 

0.70 

2130 

5.6 

63.4 

±  2.0 

0j68 

45 

11 

0.38 

5.5 

57.2 

±  5.2 

0.61 

10  July 

5.5 

50.5 

+  3.0 

0.54 

0830 

5.65 

47.1 

±6.1 

0.51 

36 

22 

0.75 

5.5 

45.1 

±  3.3 

0.48 

Exp.  Period  22-25  July  '81  Final  Ice  "thickness  hf  =  4.2 

Concentration  1.05$ 

Note:  Warm-up  started  24  July,  at  1030 

F I exura I 

Date,  Ice  "thickness  strength  (kPa)  a/0\ 

time  h(cm)  a  a'  a'/a  (0|  ■*  93.1  kPa) 


93.8  ±  5.4 

78.8  +  2.1 


Elastic  Warm- 

modulus  (MPa)  up  time 


i 

V 

W. 

24  July 
0830 

4.15 

4.10 

4.10 

82.1  + 
70  J5  + 
80.9  + 

5.5 

4.0 

2.9 

51 

24  July 

4.2  -  4.3 

59.9  + 

4.1 

0.71 

m 

1200 

4.2 

56.6  ± 

0.4 

0.67 

1.5  0.09 

ij 

4.0 

72.2  ± 

10.9 

0.85 

/;  . 

i#‘  • 

K’j 

24  July 

4.3 

30.2  + 

3.1 

0.36 

2300 

4.3  -  4.4 

21.6  + 

4.3 

0.26 

12.5  0.709 

4.25 

23.6  ± 

6.9 

0.28 

M 

25  July 

4.0 

15.1  + 

2.0 

0.18 

1000 

4.1 

13.2  + 

2.3 

0.16 

23.5  1.33 

»  * 

■* » 

4.2  -  4.3 

15.3  + 

3.3 

0.18 

*  / 

's 

Exp.  Period 

18-21  Aug 

>81 

Final  Ice 

1h  Ickness  hf  *  2.50  cm 

'~s 

Concentration  0.45$ 

Note:  Warm-up  started 

19  Aug., 

at  1530. 

* 

Flexural 

E  las  tic 

Warm- 

Date, 

Ice  Iti  Ickness  strength  (kPa) 

a/ at  modulus  (MPa) 

0 

up  time  t/h 

•  • 

time 

h  (cm) 

0  o’ 

a'/a  ( 0|  *  in  kPa)  E  E’ 

t(hr  >  f 

m 

19  Aug. 

2.15 

76.2  + 

10.1 

408 

0900  -  1000 

ft 

19  Aug. 

2.45 

76.7  + 

3.1  74.9  ±  1.1 

0.98  0.69  192 

1.5  0.24 

1630  -  1700 

19  Aug.  2.50  47.2  +6.9  32.3  +  4.6  0.68 

9 inn  -  99xn 


6.5  1.04 


Date, 

time 

lea  it  lekness 
h  (cm) 

Flexural 
strength  (kPa) 
a  o' 

o'/o 

O/Ol 

(a,  •  75.2  kPa) 

Elastic 

modu  1  us  (MPa ) 

E  E* 

Warm¬ 
up  time 
t(hr  > 

CM 

$ 

20  Aug. 

2.50 

00 

CM 

±  2.4 

16.9 

±  0.8 

0.59  0.25 

69 

17.5 

2.80 

0830  -  0930 

20  Aug. 

2.50 

19.2 

i  2.2 

17.8 

±  1.9 

0.93  0.17 

23.0 

3.68 

1430 

21  Aug. 

2.50 

24.7 

*  3.5 

NA 

(refreeze)  NA 

0500 

21  Aug. 

2.6 

36  Jo 

±  2.9 

NA 

(refreeze)  NA 

2.6 

61.7 

±  2.9 

NA 

62 

NA 

Exp.  Period  22-28  Aug  '81  Final  Ice  1h  lekness  hf  =  4.3  cm 

Concentration  0.45f 

Note:  Warm-up  started  26  Aug.,  at  1515 


Flexural 

Elastic 

Warm- 

Date, 

Ice  1h  lekness 

strength  (kPa) 

o/0| 

modulus  (MPa) 

up  time 

t/h2 

time 

h  (cm) 

a 

o' 

o'/o  (0| 

«  226  kPa ) 

E  E' 

t(hr) 

f 

25  Aug. 

1.70 

126.5 

±  10.2 

85.2 

±  2.8 

0.6  7 

270 

2015  -  2030 

25  Aug. 

3.25 

167.1 

±  10.3 

104.9 

±  2.8 

0j63 

337 

1100 

26  Aug. 

3.8 

213.6 

±  18.4 

1510 

3.56 

177.2 

±  5.5 

81.3 

±6.4 

0.46 

340 

3.2 

220.0 

±  20.5 

144.0 

±  11.3 

0j6  5 

26  Aug. 

112.4 

±  9.9 

0.55 

2310 

98.6 

±  5.54 

62.0 

±  5.73 

0.6  3 

0.48 

185 

8 

0.43 

27  Aug. 

4.3 

69.95  ±  5.3 

42.5 

±  12.8 

0^1 

0845 

4.20 

6  1 J6 

±  4  j6 

27.7 

±  2.8 

0.45 

27  Aug. 

HR 

is 

mm 

1610 

Rm i 

Em 

e xsm 

IKS 

64 

25 

1.35 

27  Aug. 

4.3 

46.8 

±  4.7 

27.7 

±  5.4 

0.59 

0.23 

2200 

4.35 

46.7 

±  5.4 

16  .7 

±  2.7 

0.36 

0.23 

53 

31 

1.68 

37 


Exp.  Period  29  Auq  -  2  Sept  '81  Final  Ice  thickness  ht  *  5.7  cm 

Concentration  0.45f 

Note:  Warm-up  started  31  Auq.,  at  1000 


a 

►  I 

a 


Ice  Iti  I ckn ess 
h  (cm) 


Flexural 
strength  (kPa) 


Elastic 
modu I  us  (MPa ) 
kPa)  E  E ' 


up  time 
t(hr ) 


30  Aug.  2.9  181.0  ±  8.8  79.8  +  7.0 

0900  -  1005  2.80  174.0  ±  7.4  66.9  ±  7.8 

2.55  -  2.60  165.0  ±  12.1  82j6  ±  10.8 


31  Aug. 

0810  -  0930 


31  Auq. 

5.5 

1 

1630  -  1730 

540 

5.5 

167.5  +  13.6 

170.5  ±  8.9 
163.2  i  1.8 


107.1  ±  1.7  0.64 

103.1  +  8.9  040 

89.1  ±  4.7  0.55 


1  Sept. 
2230 


2  Sept. 
0840 


.23 

.23 

49 

37.5 

1.12 

Exp.  Period  3-9  Sept  '81  Final 

Concentration  0.70J 

Note:  Warm-up  started  7  Sept.,  at  0900 

Flexural 


Final  Ice  Thickness  hj  •  6,5  cm 


Ice  thickness 
h  (cm) 


strenqth  (kPa) 
0  O' 


6  Sept. 

1730  -  1830 


3.83  -  3.90  148.6  t  4.0  85.2  ±  3.4 

3.8  -  3.9  124.7  ±  12.3  77.7  ♦ 

128.0  +  5.1  108.1  ±6 


O'/O  (0| 


Elastic 

/01  modulus  (MPa) 

\76  kPa )  E  E' 


up  time  t/h 
t(hr)  t 


1.52 

146 

.« 

187 

176 

1.57 

142 

169 

192 

1.85 

7  Sept. 

0830  -  0850  54 


151.9  +  5.5 
191.4  ±  2.2 


80.3  t  5.5 


179  182 


67.75  +  1.4  043 

61.0  t  0.1  040 

624  ±  1.0  046 


95 


99 


I  cm  thickness 
h  (cm) 


flexural 
strength  <kPa) 


8  Sept. 

0800  -  0900 


75.6  +  1.7 

80.3  +  1j5 

75.8  +  1.4 


6.5  44.7  +  5.6 
6.50  -  6.55  40.5  +  6.0 
6.15  -  6.25  44.8  +  5.9 


33.4  ±  2J6 
36  .0  +  2.2 
39.1  +  1.2 


18.9  +  0.9 
20.5  +  2.8 
17.4  +  1.7 


Elastic  Warm- 

2 

0/0 j  modulus  (MPa)  uo  time  t/h 

o'/o  (0|  •  175  kPa)  E  E'  t(hr)  f 


43  23.5  0.62 


Exp.  Period  10  -  14  Sept  '81  Final  ice  Ihlekness  hf  *  3.6  cm 

Concentration  0.70$ 

Note:  Warm-up  started  11  Sept.,  at  1000 

Flexural 

Date,  Ice  Ihlekness  strength  (kPa)  o 

time  h  (cm)  O  O'  O'/O  <0|  * 


10  Sept 


11  Sept. 
1200  -  1315 


11  Sept.  2. 

1630  -  1715  3. 

3.45 


125.1  ±  9.0 

118.2  +  9.5 
130.7  +  3.4 


105.1  +  3.0 

91.1  +  4.9 
100,6  +  5J5 


75.6  +  3.5 
5  7.4  +  9.7 
72.5  +  2.5 


62.2  +  3.1 
6  5,6  +6.7 
65.8  +  2.7 


63.9  +  2.2 
45.7  +  3.3 

62.9  +  2.1 


36.4  +  5.8 
32.8  +  3.2 

28.4  +  2.4 


Elastic  Warm- 

2 

o/oi  modulus  (MPa)  up  time  t/h 

O’/O  (0|  «  143  kPa >  E  E*  t(hr)  f 


3 

4 

111 

124 

2 

0.16 

5E 

0 

12  Sept. 
0100 


32.0  ±  1£ 


1 1  j6  +  0.7  0.35 


15  1.16 


12  Sept.  3.55 

0900  -  1020  3j5 

3.55 


40.4  +  4.3 
29.0  +  2.0 
25.7+2.8 


20.3  +  1.8 
16.9  +  1.3 
15.6  +  1.8 


35  37  23.5  1.81 


12  Sept. 
1900  -  2045 


28.3  t  3.0 
30.8  +  1  J5 

27.3  ±  1.5 


10.1  +  0.1 

13.2  +  0.7 
18.5  +  1.1 


22  23  34  2j62 


13  Sept. 
1740 


17.2  +  0.3 


13.9  +  2J5  0.81 


55.5  4.28 


Exp.  Period  14-18  Sept  '81  Final  le 

Concentration  0.70 t 

Note:  Warm-up  started  15  Sept.,  at  0830 

Flexural 

Date,  Ice  1h  Ickness  strength  (kPa) 

time  n  (cm)  a  o' 


Final  Ice  thickness  hf  *  2.3  cm 


15  Sept. 

2.10 

81.3 

+  8.7 

46.0  +  5.6 

0.57  ( 

0900  -  0950 

2.05 

88,8 

*  2.8 

47.2  +  5.2 

0.53  ( 

• 

2.00 

91.2 

+  9.8 

57.9  +  0.7 

0  J54  ( 

j  15  Sept. 

2.2 

59.6 

+  14.7 

47.0  t  3.0 

0.79  1 

1  1120  -  1215 

2.2 

62.8 

*  2.4 

44.3  ±  3.3 

0.70  ( 

2.1 

64.4 

+  2.4 

40.1  +  4.4 

0j52  ( 

Elastic  Warm- 

a/ a  I  modulus  (MPa)  up  time  t/h 

a'/0  (0|  »  114  kPa)  E  E-  t(hr)  t 


.71 

.78 

.80 

69 

0.5 

0.09 

15  Sept. 
1700  -  1730 


37.3  +  2.6 
38.0  +  7.5 


25.7  0.6  9 

28.7  +  2.1  0.* 
25.4  ♦  5.7  0j68 


16  Sept. 

2.2  -  2.3 

17.9 

+  0.7 

16.3  +  1.1 

0.91  I 

0900  -  1010 

2.2 

14.2 

*  1/5 

14.7  *  2.4 

1.04  1 

2.1 

14.8 

+  2.1 

( 

Exp.  Period  19  -  21  Sept  '81  Final  Ic. 

Concentration  0.70 1 

Note:  Warm-up  started  20  Sept.,  at  0700 

Flexural 

Date,  Ice  thickness  strenqlh  <kPa> 

time  h(cm)  a  O’ 


20  Sept. 
0745  -  0850 


Final  Ice  it  Ickness  ht  *  3.4  cm 


0/  Oi  imiu 

(j Wo  (<J|  *  135  kPo)  E 


modulus  (MPa)  up  time  t/h' 


1  61.8 

±  3.6 

49.7  +  4.0 

0.81 

1  69J5 

i  3.6 

50.1  t  26 

0.62 

15  63.9 

+  3.9 

47.7  t  6 ,4 

0.75 

21  Sept. 
1430  -  1440 


3.25  16.6  +0.7 


31.5  2.72 


Exp.  Period  26  -  28  Sept  *81  Final  Ic 

Concentration  0.951 

Note:  Warm-up  started  27  Sept.,  at  0930 

F I exura I 

Date,  Ice  thickness  strenqlh  (kPa) 

time  h(cm)  <j  <j' 


Final  Ice  thickness  •  3.1 


Elastic 


a/a |  modu  I  us  (MPa  > 

O’ /a  (<J|  *  66.2  kPa)  E  E' 


Warm¬ 
up  time  t/h' 
t(hr ) 


27  Sept 


27  Sept. 
1315  -  1420 


27  Sept. 
1800  -  1835 


27  Sept. 
2200 


28  Sept. 
0910  -  0920 


67.5  +  3.6 
62.4  ±  2.0 
69.7  +  5.9 


39.7  +11.4 
43.5  ±  5.6 


21.8  +  2.7 

23.31  ±  2.3 


32.3  +  4.1  0.48 
36.1  ±  4.4  0.58 
36.9  +  4.7  0.53 


25.06  +  3.2  0.66 

24.7  ±  2.9  0.57 


13.7  +  0.0  0.6  3 
14.5  t  1.5  0.62 


23.5  ±  1.0  10.2  +  1.7  0.43 


3.05  24.4  +  5.2  25,3  +  1.2  0.74 


12.5  1.30 


Exp.  Period  16  -  21  Oct  '81  Final  Ice  It 

Concentration  0.951 

Note:  Warm-up  started  19  Oct.,  at  0930 

Flexural 

Date,  Ice  thickness  strength  (kPa) 


Final  Ice  Ihlckness  hf  *  5.0  cm 


O/Oi 

a'/a  (0|  *  130  kPa) 


Elastic  Warm- 

modulus  (MPa)  up  time  t/h* 


\  17  Oct. 

2.00 

86.9  +  7.0 

47.2  +6.4 

0.54 

im  2110  -  2135 

2.00 

91.5  +  21.5 

45.5  +  5.1 

0.50 

118 

| 

2.00 

66.5  +  3.2 

38.5  +  3.1 

0.58 

18  Oct. 

3.22  -  3.29 

114.95  +  8.3 

60.4  +  4.3 

0.52 

1340  -  1430 

3.24  -  3.27 

99.2  ±  3 j6 

56.4  +  3.1 

0.57 

131 

*, 

3.31  -  3.37 

106  J5  ±  4.4 

43.7  +  2.2 

0.41 

L  18  Oct. 

3.68  -  3.80 

121.54  +  11. 

56.1  +  4.5 

0.46 

122 

1  1950 

19  oct. 

4.58  -  4j64 

143.5  +  3.3 

59.4  +  1.8 

0.41 

158 

V-  0850 

,  e 

F, 

20  Oct. 

4.95  -  5.15 

49.2  +  3.9 

18.0  +  2.3 

0.37 

0.38 

n 

* 

1520  -  1610 

5.10 

46.1  +  2.5 

14.2  +  0.9 

0.31 

0.36 

43 

25.3  1.01 

5.05 

46.3  +  8.8 

10.5  +  0.2 

0.22 

0.36 

■ ' 

20  Oct. 

5.00 

29.8  +  3.5 

10.0  +  0.5 

0.40 

0.23 

25 

33.0  1.32 

• 

2300 

21  Oct. 

5.00 

10.8  ♦  1.8 

0.08 

3 

0020  -  0900 

5.00 

6.74  +  1.0 

0.05 

11 

42.5  1.70 

5.6  cm 


Exp.  Period  15  -  16  Nov  '81  Final  ice  "thickness  ht  = 

Concentration  0.951 


Note:  Warm-up  started 

at  1300.  Overhead 

doors 

between  prep,  tank 

and 

basin  opened 

15  Nov.  at  2300  to 

accelerate 

warm-up. 

F 1 exura 1 

E  lastic 

Warm¬ 

Date, 

1  ce  th  1  ckn  ess 

strength  (kPa) 

O/O: 

modulus  (MPa) 

up  time 

t/h  2 

time 

h  (cm) 

a  o' 

a’/a 

(0| 

=  135  kPa) 

E  E' 

t(h  r ) 

f 

1 3  Nov. 

2.1 

90.8 

0800 

2.1 

76.2 

2.3 

76.7 

14  Nov. 

5.0 

123.4 

1030 

5.1 

114.9 

5.1 

136.2 

14  Nov. 

1900 

5.4 

94.9 

0.70 

6 

0.19 

5.4 

91.2 

0.68 

1 5  Nov. 

5.4 

70  j6 

0.52 

1100 

5.6 

69.5 

0.52 

64 

22 

0.70 

5.7 

74.8 

0.55 

15  Nov. 

5.5 

71.4 

0.53 

2030 

5j6 

63.8 

0.47 

63 

31.5 

1.00 

5.75 

55.6 

0.41 

16  Nov. 

5j6 

18.0 

0.13 

0700 

5j6 

18.5 

0.14 

32 

42.0 

1.34 

5.75 

14.5 

0.11 

Exp.  Period 

18  -  19  Nov  ' 

81  Final  ice 

th  Ickness  hf  =  3.3 

cm 

Concentration  0.95< 

Note:  Warm-up  started  18 

Nov.,  at  0830 

Flexural 

E  lastic 

Warm¬ 

Date,  1 

Ice  thickness 

strength  (kPa) 

a/cTi 

modulus  (MPa) 

up  time 

t/h  2 

time 

h  (cm) 

o  o' 

a'/a 

<CT| 

»  104  kPa ) 

E  E' 

t(h  r ) 

t 

18  Nov. 

0800 

2.7 

87.3 

62 

2.8 

101.3 

18  Nov. 

2.9 

89.1 

0.86 

0.5 

18  Nov. 

3.1 

56.1 

0.54 

1200 

3.2 

86.3 

0.83 

36 

3.5 

0.32 

3.1 

74.0 

0.71 

Warm- 


3 

} 

i  • 

p 

k 

*p 

ij 

} 

! 

i 

» 

i 

l 

i 


Flexural 

Elastic 

Warm- 

Data, 

Ice  ih  Ickness 

strenqth  (kPa) 

O/Oi 

modulus  (MPa) 

up  time 

t/h  2 

time 

h  (cm) 

a  o' 

o'/o  (<7|  *  104  kPa)  £  E* 

t(h  r ) 

t 

18  Nov. 

30.5 

0.29 

2050 

0.38 

19 

12.5 

1.15 

Hifioii&iii 

BWilKliiflg 

0.31 

19  Nov. 

3.3 

23.9 

0.23 

06  50 

5.4 

22.6 

0.22 

8 

22.0 

2.02 

3.3 

27.0 

0.26 

Exp.  Period  20  -  23  Nov 

'81  Final  Ice  tfi  Ickness  hf  = 

5.6  cm 

Concentration  0.95f 

Note:  Warm-up  started  21  Nov.,  at  0900 

F 1 exura 1 

Elastic 

Warm- 

Date, 

Ice  th  ickness 

strenqth  (kPa) 

o/0\ 

modulus  (MPa) 

up  time 

t/h  2 

time 

h  (cm) 

a  o' 

a'/o  < CTj  =  135  kPa) 

E  E' 

t(h  r ) 

t 

20  Nov. 

2.4 

84.0 

1200 

2.4 

90.3 

66 

21  Nov. 

5.0 

138.8 

0830 

4.8 

149.1 

188 

4.6 

143.1 

21  Nov. 

5J6 

104.3 

0.77 

1730 

5.6 

104.9 

0.78 

66 

8.5 

0.27 

5.4 

104.5 

0.77 

22  Nov. 

5.7 

71.5 

0.53 

0930 

5.7 

61.9 

0.46 

50 

24.0 

0.77 

5.5 

54.7 

0.40 

22  Nov. 

1900 

5.7 

30.4 

0.22 

41 

34.0 

1.08 

5.5 

42.9 

0.32 

23  Nov. 

5J6 

15.9 

0.12 

0530 

5.6 

16.5 

0.12 

14 

44.5 

1.42 

5.5 

19.5 

0.14 

43 


3 


i 


Exp.  Period  28  Nov  -  1  Dec  '81  Final  Ice  16  Ickness  hf  =  5.6  cm 

concentration  0.95< 

Note:  Warm-op  started  29  Nov.,  at  1230 


F 1 exura 1 

Elastic 

Warm- 

Date, 

Ice  16  ickness 

strenqt6  (kpa) 

a/a; 

modu  1  us  (MPa ) 

up  time 

t/h  2 

time 

h  (cm) 

0  o' 

o'/a 

(0|  =  135  kPa ) 

E  E' 

t(h  r ) 

f 

28  Nov. 

2.1 

6  0.6 

46 

0900 

2.1 

6  1.0 

28  Nov. 

2.9 

98.5 

1800 

2.9 

102.4 

89 

3.0 

108.1 

29  Nov. 

4.8 

146.7 

1100 

4.8 

142.0 

4.7 

150.0 

183 

29  Nov. 

5.4 

107.1 

0.79 

1800 

5.4 

106 .1 

0.79 

135 

5.5 

0.18 

5.3 

108.3 

0.80 

30  Nov. 

5.7 

64.7 

0.48 

0830 

5.6 

58.3 

0.43 

60 

20.0 

0.64 

5.5 

61.2 

0.45 

30  Nov. 

5.6 

62.3 

0.46 

1600 

5.6 

59.6 

0.44 

60 

27.5 

0.88 

5.6 

75.7 

0.56 

30  Nov. 

5.6 

52.2 

0.39 

2200 

5.7 

55.0 

0.40 

41 

33.5 

1.07 

5.5 

52.0 

0.39 

1  Dec. 

5.6 

48.3 

0.36 

0530 

5.7 

31.5 

0.23 

33 

41.0 

1.31 

5.5 

43.9 

0.33 

1  Dec. 

0730 

5.6 

33.3 

0.25 

32 

43.0 

1.37 

5.6 

30.2 

0.22 

1 


I 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 
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Properties  of  urea-doped  ice  in  the  CRREL  test  basin  / 
by  Ken-ichi  Hirayama.  Hanover,  N.H.:  Cold  Regions  Re¬ 
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